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Abstract. The solar neutrino problem is reviewed and the possible vacuum oscillation 
and MSW solutions of the problem are considered. 



1 Introduction 

The problem of neutrino mass is the central problem of present day neutrino 
physics and one of the central problems of contemporary elementary particle 
physics. ° The existence of nonzero neutrino masses is typically correlated in the 
modern theories of the elementary particle interactions with nonconservation of 
the additive lepton charges, Le, and (see, e.g., [1]). The rather stringent 
experimental limits on neutrino masses obtained so far together with cosmolog- 
ical arguments imply (see, e.g., [7]) that if nonzero, the masses of the flavour 
neutrinos must be by many orders of magnitude smaller than the masses of the 
corresponding charged lepton and quarks belonging to the same family as the 
neutrino. The extraordinary smallness of the neutrino masses is related in the 
modern theories of electroweak interactions with massive neutrinos to the ex- 
istence of new mass scales in these theories. Thus, the studies of the neutrino 
mass problem are intimately related to the studies of the basic symmetries of 
electroweak interactions; they are also closely connected with the investigations 
of the possibility of existence of new scales in elementary particle physics. Cor- 
respondingly, the experiments searching for effects of nonzero neutrino masses 
and lepton mixing are actually testing the fundamental symmetries of the elec- 
troweak interactions. These experiments are also searching for indirect evidences 
for existence of new scales in physics. 

Neutrinos are massless particles in the standard (Glashow-Salam- Weinberg) 
theory (ST) of the electroweak interactions. The observation of effects of nonzero 



The neutrino mass problem, the phenomenological implications of the nonzero neu- 
trino mass and lepton mixing hypothesis, the properties of massive Dirac and mas- 
sive Majorana neutrinos, the neutrino mass generation in the contemporary gauge 
theories of the electroweak interaction as well the role massive neutrinos can play 
in astrophysics and cosmology are the subject of a number of review articles and 
books: see, e.g., [l]-[6]. 
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neutrino masses and lepton mixing would be a very strong indication for exis- 
tence of new physics beyond that predicted by the standard theory. The studies 
of the neutrino mass problem can lead to a progress in our understanding of the 
nature of the dark matter in the Universe as well [8] . 

One of the most interesting and beautiful phenomenological consequences 
of the nonzero neutrino mass and lepton mixing hypothesis are the oscillations 
of neutrinos [9], i.e., transitions in flight between different types of neutrinos, 
Ve ^ Vy. and/or Ve, <^ v-r and/or Vfj, ^ Vr, and antineutrinos, Pg I'm and/or 
j/e ^ v'r and/or ij^ ^ Vr- If, for example, a beam of Ve neutrinos is produced 
by some source, at certain distance R from the source the beam will acquire a 
substantial component if Ve Vfj. oscillations take place. The probability to 
find at distance R from the source of when the neutrinos propagate in 
vacuum and the massive neutrinos are relativistic, P{ve i^ti), is a function of 
the neutrino energy E, the differences of the squares of the masses of the 
neutrinos having definite mass in vacuum, zim^fe = m| — ml, and of the 
elements of the lepton mixing matrix U (sec, e.g., [1], [2]). 

At present we have several indications that neutrinos indeed take part in 
oscillations, which suggest that neutrinos have nonzero masses and that lepton 
mixing exists. One of the indications comes from the results of the LSND neutrino 
oscillation experiment performed at the Los Alamos meson factory [10]. The 
events observed in this experiment can be interpreted as being due to i?^ <-> 
Pe oscillations with Am^ ^ few eV"^ and sin^ 2^ ^ few x 10~^, where Arrv^ 
and sin^ 26 are the two parameters - the neutrino mass squared difi^erence and 
the neutrino mixing angle, which characterize the oscillations in the simplest 
case. The second indication is usually referred to as the atmospheric neutrino 
problem or anomaly [11], [12] : the ratio of the //.—like and e— like events produced 
respectively by the fluxes of (t'^ + i/^) and {ve + i^e) atmospheric neutrinos with 
energies ~ (0.2 — 10.0) GeV, detected in Kamiokande, 1MB, Soudan and Super- 
Kamiokande experiments, is smaller than the theoretically predicted ratio. The 
atmospheric neutrino data can be explained by Vfj, Vr and z?^ <-> 9t oscillations 
with Am^ ~ (10""^ — 10~^) eV"^ and a relatively large value of sin^ 29 close to 1. 

The amount of the solar neutrino data available at present, the numerous 
nontrivial checks of the functioning of the solar neutrino detectors that have 
been and are being performed, together with recent results in the field of solar 
modeling associated, in particular, with the publication of new more precise 
helioseismological data and their interpretation, suggest, however, that the most 
substantial evidence for existence of nonzero neutrino masses and lepton mixing 
comes at present from the results of the solar neutrino experiments. In view of 
this we will devote the present lectures to the solar neutrino problem and its 
possible neutrino oscillation solutions. 

The "story" of solar neutrinos begins, to our knowledge, in 1946 with the 
well-known article by B. Pontecorvo [13], published only as a report of the Chalk 
River Laboratory (in Canada). In it Pontecorvo suggested that reactors and the 
Sun are copious sources of neutrinos. On the basis of neutrino flux and inter- 
action cross-section estimates he concluded [13] that the experimental detection 
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of neutrinos emitted by a reactor (i.e., the observation of a reaction caused by 
neutrinos) is feasible, while the detection of solar neutrinos can be very difficult 
(but not impossible). In the same article the radiochemical method of detection 
of neutrinos was proposed. As a possible concrete realization of the method, 
a detector based on the Cl-Ar reaction i/^ + ^^Cl '^^Ar + e~ was discussed. 
The possibility to use the Cl-Ar method for detection of neutrinos was further 
studied in 1949 by Alvarez [14]. A Cl-Ar detector for observation of solar neutri- 
nos was eventually built by Davis and his collaborators [15]. The epic Homestake 
experiment of Davis and collaborators, in which for the first time neutrinos emit- 
ted by the Sun were detected, began to operate in 1967 and still continues to 
provide data. It was realized in 1967 as well [16] that the measurements of the 
solar neutrino flux can give unique information not only about the physical con- 
ditions and the nuc;lear reactions taking place in the central part of the Sun, but 
also about the neutrino intrinsic properties. 

The solar neutrino problem emerged in the 70'ies as a discrepancy between 
the results of the Davis et al. experiment [15], [17] and the theoretical predictions 
for the signal in this experiment [18], based on detailed solar model calculations. 
The hypothesis of unconventional behaviour of the solar on their way to the 
Earth (as like, e.g., vacuum oscillations [9], [16] '^ij,{t) and/or Uf. Vg, Vg 

being a sterile neutrino, etc.) provided a natural explanation of the deficiency of 
solar neutrinos reported by Davis et al. However, as the fraction of the solar i/g 
flux to which the experiment of Davis et al. is sensitive (neutrinos with energy E 
> 0.814 MeV) was known [18] i) to be produced in a chain of nuclear reactions 
(representing a branch of the pp cycle) which play a minor role in the physics of 
the Sun and whose cross-sections cannot all be measured directly in the relevant 
energy range on Earth, and ii) to be extremely sensitive to the predicted value 
of the central temperature, Tc, in the Sun (scaling as T^^), the possibility of an 
alternative (astrophysics, nuclear physics) explanation of the Davis et al. results 
could not be excluded. 

In 1986 an independent measurement of the high energy part (E > 7.5 MeV) 
of the flux of solar neutrinos was successfully undertaken by the Kamiokande II 
collaboration using a completely different experimental technique; in 1990 the 
measurements were continued by the Kamiokande III group with an improved 
version of the Kamiokande II detector [19]. At the beginning of the 90'ies two 
new experiments, SAGE [20] and GALLEX [21]. sensitive to the low energy part 
(E > 0.233 MeV) of the solar neutrino flux, began to operate and to provide 
qualitatively new data. The Kamiokande III detector was succeeded by an ap- 
proximately 30 times bigger version called appropriately "Super-Kamiokande" , 
which began solar and atmospheric neutrino detection on April 1, 1996 [22]. The 
data obtained since 1986 did not alleviate the solar neutrino problem - on the 
contrary, they made the case for existence of solar neutrino deficit even stronger. 

At the same time considerable efforts were also made to understand better 
the potential sources and the possible magnitude of the uncertainties in the 
theoretical predictions for the signals in the indicated solar neutrino detectors, 
and to develop improved, physically more precise solar models on the basis of 
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which the predictions are obtained [23]. Remarkable progress in this direction 
was made in the last several years with the development of the solar models 
which include the diffusion of helium and the heavy elements in the Sun [24]- 
[32] , as well as with the appearance of new more precise helioseismological data 
permitting new critical tests of the solar models to be performed [33]-[37]. 

With the accumulation of more data and the developments in the theory 
certain aspects of the solar neutrino problem changed and new aspects appeared. 

In the present lectures we shall review the current status of the solar neutrino 
problem. We shall also review the status of the neutrino physics solutions of the 
problem based on the hypotheses of vacuum oscillations [9] or of matter-enhanced 
transitions [38], [39] of solar neutrinos. 

2 The Data and the Solar Model Predictions 

We begin with a brief summary of relevant solar model predictions and of the 
solar neutrino data. According to the existing models of the Sun [23] , the solar 
flux consists of several components, six of which are relevant to our discussion: 

i) the least energetic pp neutrinos (E < 0.420 MeV, average energy E = 
0.265 MeV), 

ii) the intermediate energy monoenergetic ^Be neutrinos (E=0.862 McV (89.7% 
of the flux), 0.384 MeV (10.3% of the flux)), 

iii) the higher energy *B neutrinos (E < 14.40 MeV, E = 6.71 MeV), and three 
additional intermediate energy components, namely, 

iv) the monoenergetic pep neutrinos (E=1.442 MeV), and the continuous spec- 
trum CNO neutrinos produced in the /3+— decays 

v) of (E < 1.199 MeV, E = 0.707 MeV), and 

vi) of (E < 1.732 MeV, E = 0.997 MeV). 

The pp, pep, ''Be and *B neutrinos are produced in a set of nuclear reactions 
shown in Fig. 1. These make part of three major cycles (the pp-cycles) of nuclear 
fusion reactions in which effectively 4 protons burn into ^He with emission of two 
positrons and two neutrinos, generating approximately 98% of the solar energy: 

4p^ *He + 2e+ -I- 2!/e. (1) 

The flrst (pp-I) cycle (or chain) begins with the p-p (or p-e~-p) fusion into 
deuterium and ends with the reaction ^He+^He ^Hc + 2p. The second (pp-II) 
and the third (pp-III) cycles begin with the production of ^Be in ^He+^He fusion 
and end respectively with the processes ''Li+p 2^He and ^B — > 2^He+c+ +/^e 
(see Fig. 1). 

The CNO neutrinos are produced in the CNO-cycle of reactions, which, ac- 
cording to the present day understanding, plays minor role in the energetics of 
the Sun: + p ^ i-'^N + 7, ^-''N ^ + e+ + v,, ^'''C + p ^ + 7, 
i4n + p ^ 150 + 7, ^ + e+ + i^e and + p ^ +4 jje. The pp, 
pep, '^Be and *B neutrino spectra are depicted in Fig. 2. Let us note that the 
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(%) (MeV) 

p + p + e+ + (99.96) < 0.420 



or 



p + e- + p 2jj _^ ^0.44) 1.442 

2H + p ^ 3jje + ^ (100) 
^He + 3He ^ a + 2 p (85) 



or 
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He + *He ^ ^Be + 7 (15) 

7r . - 7t- ^ /1.x / 0-861 90% 

+ ^ + 1 0.383 10% 

■^Li + p -> 2 a 

or 

■^Be + p ^ *B + 7 (0.02) 

^B ^ ^Bg* + e + + z^, < 15 

^Be* ^2 a 
or 

^He + p ^ *He + e+ + (0.000004) 18.8 



Fig. 1. The nuclear reactions of the pp-chain in the Sun. 



shapes of the continuous spectra of the pp, B and the CNO neutrinos are to 
a high degree of accuracy solar physics independent. The total fluxes of the pp, 
pep, ''Be, *B and CNO neutrinos depend, although to a different degree, on the 
physical conditions in the central part of the Sun [23] (see further) . 

Solar neutrinos are produced in the central solar region (which practically 
coincides with energy production region) with radius Th = 0.25 Rq, Rq = 6.96 x 
10^ km being the radius of the Sun. The dependence of the source-strength 
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Fig. 2. The spectra of the pp. pop, ^Bc and **B neutrino fluxes (from [23]). Shown are 
also the ranges of solar neutrino energies to which the Ga - Ge, CI - Ar and Kamiokande 
II and III experiments are sensitive. 



functions for the pp, ^Be and neutrinos on the distance from the center 
of the Sun, r, is shown in Fig. 3. As this figure illustrates, the major part of 
the *B neutrinos flitx is generated in a rather small region, r ^ 0.10 Rq close 
to the center of the Sun; the region of production of ''Be neutrinos extends 
to r = 0.15 Rq, while the region of the pp neutrino production is the largest 
extending to r = 0.25 Rq. 

Three different methods of solar neutrino detection have been and are being 
used in the six solar neutrino experiments [17], [19]-[22] that have provided data 
so far: the Cl-Ar method proposed by Pontecorvo [13] - in the experiment of 
Davis et al. [15], [17], the u — e~ elastic scattering reaction in the (Super-) 
Kamiokande experiments [19], [22], and the radiochemical Ga-Ge method - in 
SAGE [20] and GALLEX [21] experiments. 

The threshold energy of the reaction I'e CI e~ +^'^ Ar on which the 
Cl-Ar method is based, is Eth(Cl) = 0.814 MeV. Consequently, the pp neutrinos 
do not give contributions in the signal in the Davis et al. detector. Inspecting 
the predictions of all solar models presently discussed in the literature one finds 
that the major contribution to the signal in the Cl-Ar experiment, between 64% 
and 79%, should be due to the ^B neutrinos; the ^Be neutrinos are predicted to 
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Fig. 3. The production of tiie fluxes of pp, ^Be and *B neutrinos as a function of tlie 
distance R from the center of the Sun (the latter is expressed in units of the solar 
radius R©) (from [23]). 



generate between 22% and 13% of the total signal, and the pep and the CNO 
neutrinos - between 13% and 8%. With a threshold neutrino energy, Eth(K), 
first of 9.5 MeV (7.5 MeV) and subsequently reduced to 7.5 MeV and further to 
7.3 MeV (6.5 MeV), the (Super-) Kamiokande experiments can detect only the 
higher energy component of the solar v,. flux. 

Having the lowest threshold energy Eth(Ga) — 0.233 MeV, the Ga-Ge detec- 
tors GALLEX and SAGE are sensitive to all six components of the flux consid- 
ered above. Moreover, the major part of the signal in these detectors, between 
51% and 63%, is predicted to be produced by the pp neutrinos; the ^Be neutri- 
nos are expected to generate between 28% and 24% of the total signal, the ^B 
neutrinos - between 12% and 5%, and the CNO neutrinos ~ around (5~7)%. 

The above analysis implies that the Cl-Ar and Kamiokande experiments on 
one side, and the Ga-Ge experiments on the other, are most sensitive to very 
different components of the solar neutrino flux: the former - to the ^B neutrinos, 
and the latter - to the pp neutrinos. The ^Bc neutrinos are predicted to give the 
second largest (and non-negligible) contributions to the signals in both Cl-Ar 
and Ga-Ge experiments. 

Let us turn next to the data. The average rate of ^^Ar production by solar 
neutrinos, R(Ar), observed in the experiment of Davis et al. in the period 1971- 
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1996 (altogether 800 solar induced events registered) is [17] 

(2.56 ±0.16 ±0.14) SNU. (1) 

Here (and in the experimental results we quote further) the first error is statis- 
tical (1 s. d.) and the second error is systematic. The flux of neutrinos, ^-q, 
measured by the Kamiokandc experiments reads [19] 

= (2.80 ± 0.19 ± 0.33) X 10^ cm-^sec'^ (2) 

The result is based on a statistics of about 600 events accumulated by the two 
experiments in 2079 days of measurements in the period 1986 - 1996. 

The GALLEX and SAGE experiments began to collect data in 1991 and 
1990, respectively. The GALLEX group has registered (in 65 runs) approximately 
300, and the SAGE group has registered (in 33 runs) about 100 solar neutrino 
induced events. The average rates of ^^Ge production by solar neutrinos, R(Ge), 
measured by the two collaborations are [21], [20] 

RGALLEx(Ge) = (76.2 ± 6.5 ± 5) SNU, (3) 

RsAGE(Ge) = (73 ± 8.5 1^:^) SNU. (4) 

Obviously, the results of the two experiments are compatible. Adding the statis- 
tical and the systematical errors in (3) and in (4) in quadratures and combining 
the two results (i.e., taking the weighted average) we find 

Rcxp(Gc) = (75.4 ± 7) SNU. (5) 

Recently, the Super-Kamiokande collaboration announced results on solar 
neutrinos based on data collected during a period of 306.3 days. About 4000 (!) 
solar neutrino induced events have been detected. The following value of the 
neutrino flux was measured [22]: 

= (2.44 ± 0.06 tlf^) X 10^ cm"2sec"\ (6) 

The GALLEX collaboration has successfully performed in 1994 and in 1997 
very important (and rather spectacular) calibration experiments with an artifi- 
cially prepared powerful ^^Cr source of monoenergetic v^. (four lines: E = 746 
keV (81%), 751 kcV (9%), 426 keV (9%) and 431 keV (1%)) of known intensity 
[21]. At the beginning of the two exposures the signal due to the ^^Cr neutrinos 
was approximately 15 times bigger than the signal due to the solar neutrinos. 
The results of the experiments showed, in particular, that the efficiency of extrac- 
tion of the ''^Ge, produced by neutrinos, from the tank of the detector coincides 
(within the 10% error) with the calculated one. 

Similar calibration experiment has been successfully completed in 1996 also 
by the SAGE collaboration. These experiments demonstrated, in particular, that 
the Ga-Ge detectors are capable of detecting the intermediate energy ^Bc neu- 
trinos with a high efficiency. They represent a solid proof that the data on the 
solar neutrinos provided by the two detectors are correct. They also represent 
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the first real proof of the feasibihty of the radiochemical method invented by 
Pontecorvo [13] for detection and quantitative study of solar neutrinos. 

An extensive program of calibration studies of the Super-Kamiokande de- 
tector is presently being completed. The aim is to reach an accuracy of 1% in 
the measurement of the energy of the recoil e~ from the solar neutrino induced 
reaction v + e~ ^ v + e~ . 

3 The Data Versus the Solar Model Predictions 
3.1 Modeling the Sun 

The results of the solar neutrino experiments have to be compared with the 
corresponding theoretical predictions. Many authors have worked (and many 
continue to work) in the field of solar modeling and have produced predictions 
for the values of the pp, ''Be, ^B, pep and CNO neutrino fluxes, and for the 
signals in the solar neutrino detectors: a rather detailed review of the results 
obtained by different authors prior 1992 and the corresponding references can be 
found in [24]. The articles [25], [27]-[32] describe some of the models proposed 
after 1992. Most persistently solar models with increasing sophistication and 
precision, aiming to account for and/or reproduce with sufficient accuracy the 
physical conditions and the possible processes taking place in the inner parts 
of the Sun have been developed starting from 1964 by John Bahcall and his 
collaborators [23]. 

The solar models are based on the standard assumptions of hydrostatic equi- 
librium and energy conservation made in the theory of stellar evolution. Several 
additional ingredients are needed to determine the physical structure of the Sun 
and its evolution in time [23], [31]: 

1. the initial chemical composition, 

2. the equation of state, 

3. the rate of energy production per unit mass as a function of the density p, 
temperature T and chemical composition fx, 

4. the radiative opacity k as a function of the same three quantities, and 

5. the mechanism of energy transport. 

The initial chemical composition of the Sun is, of course, unknown. However, 
the relative abundances of the heavy elements in the initial Sun, with the excep- 
tion of the noble gases and C, N, and O, are expected to be approximately equal 
to those found in the type I carbonaceous chondrite meteorites (see, e.g., [23]). 
Using the meteoritic abundances and the measured abundances in the solar pho- 
tosphere which includes hydrogen, and taking into account the possible change 
of the heavy element abundances during the evolution of the Sun, permits to fix 
the present day ratio of the heavy element (Z) and hydrogen (X) mass fractions, 
Z/X. The knowledge of Z/X, the normalization condition X + Y + Z = 1, where 
Y is the helium mass fraction, and the requirement that the solar model repro- 
duces correctly the measured value of the solar luminosity (see further) allows to 
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determine the absolute values of the initial solar element abundances (for further 
details see [23]). 

The equation of state requires the knowledge of the degree of ionization and 
the population of the excited states of all elements present in the Sun. Stellar 
plasma effects which introduce deviations from the perfect gas law have to be 
taken into account as well. As we have mentioned earlier, the energy is generated 
in the Sun in four cycles (or chains) of nuclear fusion reactions in which effectively 
4 protons burn into ^He: 4p — > ^He + 2e+ + 2i/e. Collective plasma and screening 
effects have to be accounted for in the calculation of the corresponding nuclear 
reaction cross-sections. 

The radiative opacity k is determined by the photon mean free path. It 
controls the temperature gradient (and therefore the energy flow) in the radiative 
zone. The calculation of k requires a detailed knowledge of all atomic levels in the 
solar interior as well the cross-sections of photon scattering (elastic and inelastic), 
emission, absorption, inverse bremsstrahlung, etc. and is a rather complicated 
task. The energy transport is assumed to be by radiation in the inner part of 
the Sun, and by convection in the outer region. The border region between the 
radiative and convective zones is located at r ~ 0.7 Rq, where r is the distance 
from the solar center. 

A solar model should reproduce the observed physical characteristics of the 
Sun: the mass [42] M© = (1.98892 ± 0.00025) x 10^^ g, the present radius 
Rq = (6.9596 ± 0.0007) x lO^^ km and luminosity [29] (see also [43]) Lq = 
3.844(1 ± 0.004) X 10^^ erg s~^, as well as the measured relative photospheric 
mass abundances of the elements heavier than ^He, Z, and of the hydrogen, X: 
{Z/X)photo = 0.0245 (1 ± 0.061); actually, these quantities are used as input in 
the relevant computer calculations. An important constraint is the age of the 
Sun: T© = (4.57 ± 0.01) x 10^ yr. In order to develop a solar model one typically 
studies the evolution of an initially homogeneous Sun, having a mass Mq dur- 
ing a period of time r©. To reproduce the values of Rq, Lq and {Z j X)^hoto at 
time t = tq three parameters in the calculations are used: the initial helium and 
heavy element abundances Y and Z, and a parameter characterizing the convec- 
tion efficiency in the outer region of the Sun (the mixing length parameter). The 
latter is constrained by the value of Rq . It is usually assumed that the Sun is 
spherically symmetric. 

It should be clear from the above brief discussion that the solar modeling 
requires a rather good knowledge of several branches of physics: astrophysics, 
atomic, nuclear (elementary particle) and plasma physics. The most recent and 
sophisticated standard solar models [28]-[32], [35] include the effects of the slow 
diffusion (relative to hydrogen) of helium and the heavier elements from the 
surface towards the center of the Sun (caused by the stronger gravitational pull 
of these elements relative to hydrogen). 
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3.2 Helioseismological Constraints on Solar Models 

At present one of the most stringent constraints on the solar models are obtained 
from the helioseismological data. It was discovered experimentally as early as in 
1960 [33] that the surface of the Sun is oscillating with periods which vary in the 
interval between about 15 and 3 minutes (these oscillations are usually referred 
to as the "5 minute oscillations"). In later studies about 10^ individual oscillation 
modes have been identified experimentally and their frequencies were measured 
with an accuracy of 1 part in 10^ or better. 

The Sun's surface oscillations reflect the existence of standing pressure waves 
(p- waves) in the interior of the Sun (see, e.g., [33]). Some of these waves penetrate 
deep into the region of neutrino production. The p-mode frequencies depend on 
the physical conditions in the interior of the Sun. Using a specially developed 
inversion technique and the more precise helioseismological data on the low- 
frequency oscillations which became available recently, it was possible to recon- 
struct with a remarkable accuracy the sound speed distribution, c(r), in a large 
region of the Sun [33], [35], [36] extending from r ^ 0.05 Rq to r ^ 0.95 Rq. 
Using the same data permitted to determine the location of the bottom of the 
convective zone, ri,, and the matter density at the bottom of the convective zone, 
pb, as well [37]: 

n = (0.708 - 0.714) Rq , (7) 

P6 = (0.185 - 0.199) g/cm^ . (8) 

The implications of the helioseismological data for the solar modeling are il- 
lustrated in Fig. 4 (taken from [36]), where the ratio (csMir) — CHs{r)) / CHs{r), 
CHsif) ^'^d cs Aiir) being the sound speed distributions extracted from the helio- 
seismological data and predicted by a given solar model, is plotted for two solar 
models - without and including heavy element diffusion [29] . Only the statistical 
errors in the determination of chs (t") are shown (they are so small that they are 
barely seen) , but the general conclusions which can be inferred from such a com- 
parison remain valid after the inclusion of conservatively estimated systematical 
errors [37]. As Fig. 4 illustrates, the difference between CHsir) and csAiir) for 
the model without heavy element diffusion is so large that this model is prac- 
tically ruled out by the helioseismological data. Actually, the same conclusion 
is valid for all existing solar models without heavy element diffusion (e.g., the 
models [24], [25], [26]). 

Further studies [37] have shown, in particular, that models which have been 
especially designed to explain the observed deficiency of neutrinos by lowering 
the temperature in the central region of the Sun (see further), i.e., the so-called 
models with "mixed solar core" , also do not pass the helioseismological data test. 
Thus, of the large number of solar models proposed so far only the models which 
include diffusion of the heavy elements are compatible with the helioseismological 
data. 

The agreement of the predictions of the models with heavy element diffusion 
for c(r) with the sound speed distribution deduced from the data is quite impres- 
sive. As Fig. 4 indicates, the root mean square (r.m.s.) deviation from CHs{r) 
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Fig. 4. The ratio {csm — chs)/chs, chs and csm being the sound speed distributions 
extracted from the helioseismological data and predicted by a given solar model, as a 
function of the distance R from the center of the Sun (expressed in units of the solar 
radius Rq) (from [36]). The ratio is plotted for two solar models - without (dashed 
line) and including (solid line) heavy element diffusion [29]. Only the statistical errors 
in the determination of CHsij') are included in the analysis (they are too small to be 
seen) . 



for the model of [29] in the entire range of CHs{f) determination is ^ 0.1%. For 
the model with heavy element diffusion of [35] the r.m.s. discrepancy is ^ 0.2%. 
Adding the estimated systematic uncertainties in the determination of CHs{r) 
allows the difference \csm{t) — CHs(r)[ to be somewhat larger [37] : the ratio 
\csM{'r) — chs{''')\I CHs{f) for the model [29], however, does not exceed approx- 
imately 0.4% in the region O.2i?0 < r < O.65i?0, and about 1% in the neutrino 
production region O.OSi?© < r < 0.20i?Q, for which the helioseismological data 
is less accurate. 

As can be shown [36], in the interior of the Sun one has: c?{r) ^ T(r)/ fj,(r), 
where T is the temperature and fi is the mean molecular weight. Consequently, 
even small deviations of the solar model predictions for T and /i from their 
actual values in the Sun, 6T and Sfj,, would lead to a relatively large discrepancy 
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between the predicted an measured values of c(r): bcjc = 0.5{5T/T — Sji/ ji). 
Since, according by the standard solar models, T and ^ vary by very different 
factors in the energy (neutrino) production region,'' namely, by a factor of 1.9 
and by 39%, it is quite unlikely that the discrepancies between the predicted and 
actual values of T and ji mutually cancel to produce the remarkable agreement 
for the models with heavy element diffusion, illustrated in Fig. 4. Barring such 
a cancellation, the comparison between chs{t) and csM{f) suggests that for 
the indicated models \5T\/T, ^ 2% in the interior of the Sun, and is 

considerably smaller in most part of the Sun. Since the solar neutrino fluxes 
scale approximately as T" with n = —1.1; 8; 24 respectively for the pp, ^Bc and 
*B neutrinos [40] , Tc being the central temperature predicted by the model, the 
above results implies, in particular, that it is impossible to reduce considerably 
even the *^B neutrino flux by changing the central temperature within the limits 
following from the helioseismological data. 

3.3 Predictions for the Neutrino Fluxes and the Signals in the 
Solar Neutrino Detectors 

We shall present here the results obtained in four models [29]-[32] with heavy 
element diffusion, which can be characterized by their predictions for the total 
flux of ^B neutrinos as relatively "high flux" [29], [30], "intermediate flux" [31] 
and "low flux" [32] models. The predictions for the ^B neutrino flux and for the 
signals in the solar neutrino detectors in these models determine the correspond- 
ing intervals in which the results of practically all contemporary solar models 
compatible with the helioseismological data [34]-[37] he (see, e.g., [28], [35]). 
Thus, they give an idea about the dispersion and the possible uncertainties in 
the predictions. 

In Table 1 we have collected the results of the models of Bahcall-Pinsonneault 
(BP'95) from 1995 [29], Richard et al. (RVCD) [30], Castellani et al. (CDFLR) 
[31], and of Dar and Shaviv from 1996 (DS'96) [32] for the values of the fluxes of 
the pp, ""Be, ^B, pep and CNO neutrinos at the Earth surface. We have included 
also the estimated 1 s.d. uncertainties in the predictions for the fluxes made 
by Bahcall and Pinsonneault for their model. In Tables 2 and 3 we give the 
predictions for the contributions of each of the indicated six fluxes to the signals 
in the Cl-Ar [17] and the Ga Gc [20], [21] experiments, respectively, and quote 
the predictions for the total signals in these experiments (including the estimated 
1 s.d. uncertainty in the predictions whenever it is given by the authors). 

A comparison between the experimental results (1) ~ (6) and the correspond- 
ing predictions given in Tables 2 and 3 leads to the conclusion that none of the 
solar models proposed so far provides a satisfactory description of the solar neu- 
trino data: the predictions typically exceed the observations. This is one of the 
current aspects of the solar neutrino problem. 

' They vary respectively by the factor of 53 and by 43% in the entire region of the 
c(r) helioseismological determination. 



The Solar Neutrino Problem 13 



Table 1. Solar neutrino fluxes at the Earth surface (in units of cm ^sec ^) 
predicted by the solar models [29]-[32]. 



Flux 


BP'95 


RVCD 


CDFLR 


DS'96 


#pp X 10-10 

^pcp X 10-8 
^>Bc X 10-9 

X lO-'^ 
#N X 10-8 

X 10-8 


5.91(1 tVol) 

1.40(1 tr^i) 

5.15(1 

6.62(1 tri7) 
6.18(1 ISio) 

5.45(1 tr2i) 


5.94 

1.38 
4.80 
6.33 
5.59 
4.81 


5.99 
1.40 
4.49 
5.16 
5.30 
4.50 


6.10 
1.43 
3.71 
2.49 
3.82 
3.74 


Table 2. Signals (in SNU) in Cl-Ar detectors due to the solar neutrinos, pre- 
dicted by the solar models [29]- [32]. 


Type of neutrinos BP'95 


RVCD 


CDFLR 


DS'96 


PP 
pep 
^Be 
8B 

13N 
150 


0.0 

0.22 
1.24 
7.36 
0.11 
0.37 


0.0 

0.22 
1.15 
6.71 
0.09 
0.32 


0.0 
0.22 
1.08 
5.74 
0.09 
0.31 


0.0 

0.24 
0.89 
2.64 
0.06 
0.25 


Total: 


9.3 Hi 


8.5 


7.4 


4.1 ± 1.2 



Taking into account the estimated uncertainties in the theoretical predictions 
and the experimental errors in (1) - (6) one finds that the differences between 
the predictions and the observations are largest for the "high flux" models: the 
measured values of R(Ar), ^(B) and R(Ge) in the Cl-Ar, Kamiokande and Ga- 
Ge experiments are at least by (3.5-4.0) s.d. smaller than the predicted ones 
in [29], [30]. The "low flux" model of Dar and Shaviv reproduces the result of 
the Kamiokande and Super-Kamiokande experiments for ^^(B). However, the 
prediction of the model for R(Ge) is respectively by at least 3 s.d. higher than 
the GALLEX and SAGE results (3) and (4). The discrepancy between the solar 
model predictions and the Ga-Ge data is somewhat larger if one compares the 
predictions with the combined result (5) of the GALLEX and SAGE experi- 
ments. 

Let us discuss in somewhat greater detail the results of the four representa- 
tive models [29]-[32] for the fluxes of the pp, pep, ^Be, ^B and CNO neutrinos 
shown in Table 1. The predictions for the values of the pp and pep neutrino 
fluxes are remarkably coherent: they very from model to model at most by 3% 



14 S.T. Petcov 



Table 3. Signals (in SNU) in Ga-Ge detectors due to the solar neutrinos, pre- 
dicted by the solar models [29]-[32] . 



Type of neutrinos 


BP'95 


RVCD 


CDFLR 


DS'96 


PP 


69.7 


70.1 


70.7 


72.0 


pep 


3.0 


3.0 


3.0 


3.0 


^Be 


37.7 


35.3 


32.9 


27.1 


8B 


16.1 


15.4 


12.6 


6.1 


13N 


3.8 


3.5 


3.3 


2.4 


150 


6.3 


5.6 


5.2 


4.6 


Total: 


137 tr 


133 


128 


115 ± 6.0 



and 3.5%, respectively. Actually, the two fluxes are related [23]: ^pcp is pro- 
portional to «?pp and the coefficient of proportionality is practically solar model 
independent, being determined by the ratio of the cross-sections of the reac- 
tions p-|-e~+p^D + fo and p + p ^ D + e+ -I- fe in which the pep and pp 
neutrinos are produced in the Sun. One has [23]-[32] 

^pep = (2.3 - 2.4) X 10-3^pp. (9) 

The value of the pp flux is constrained by the existing rather precise data 
on the solar luminosity, Lq. Indeed [41], the solar luminosity is determined by 
the thermal energy released in the Sun in the two well known cycles of nuclear 
reactions, the pp (Fig. 1) and the CNO cycles (sec, e.g., [23]), in which 4 protons 
are converted into ^He with emission of 2 neutrinos. From the point of view of 
the energy eff'ectively generated, the indicated hydrogen burning reactions can 
be written as 

4p + 2e- ^ %e + 2i^e- (10) 

Depending on the cycle, the two emitted neutrinos can be both of the pp or 
pep type, or a pp (pep) and a ^Be, a pp (pep) and a *B (pp cycle), and a ^^O 
and/or ^'^N (CNO cycle) neutrinos [23], [31]. The thermal energy released per 
one produced pp, pep, ^Be and ^B neutrino is equal to (Q/2 — Ej), where Q = 
26.732 MeV is the Q-value of the reaction (7), and Ej is the average energy of the 
neutrino of the type j (j = pp, ...). The energy released per one ^^O and/or one 
^^N neutrino, as can be shown (taking into account, in particular, the discussion 
of the rates of the different reactions of the CNO cycle given in [23]) is equal with 
a high precision to the difference (Q/2 — (En -I- Eo)/2). Obviously, the values of 
Q and Ej are solar physics (and therefore solar model) independent. 

Given the average energies Ej carried away by the pp, pep, ^Be, ^B and 
CNO neutrinos (they are listed at the beginning of Sect. 2, Ese = 0.813 MeV), 
and knowing that the energy emission by the Sun is quasi-stationary (steady 
state), it is possible to relate Lq with the pp, pep, ''Be, ^B and CNO neutrino 
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luminosities of the Sun. One finds in this way the following constraint on the 
solar neutrino fluxes: 

^pp+0.958^Be+0.955#cNO+0.910^pep = (6.514 ± 0.031) xlO^" cm-^sec-^ , 

(11) 

where 

we have used [29] (see also [42], [43]) Lq = (3.844 ± 0.015) x 10^3 erg sec" ^ 
and have neglected the terms proportional to ??b and to (^n — ^o) in the left 
hand side of the equation, which are predicted to be considerably smaller than 

3 X 10® cm~^sec~^. The fluxes <?b and (<?n — ^o) have to be more than 46 
and (3 - 4) times bigger than the largest ("high flux") model predictions [29] 
and [30], respectively, in order for these terms to exceed the indicated value. 
The coefficient multiplying the <?Be term in (11) is just the ratio of the thermal 
energies produced per one ''Be and one pp neutrino, (Q/2 — EBe)/(Q/2 — Epp), 
etc. Since, as Table 1 shows, <?bc and <?cno arc smaller than at least by the 
factors 0.09 and 0.02, respectively, and ^pep is even smaller (see (9)), (11) limits 
primarily the pp neutrino flux. 

Comparing the experimental results (3) (6) with the solar model predictions 
given in Table 3 one notices, in particular, that the rate of '^^Ge production due 
only to the pp neutrinos, RPP(Ge) = (70 — 72) SNU, is very close to the rates 
observed in the GALLEX and SAGE experiments. This suggests that a large 
fraction, i.e., roughly at least half, of the pp (electron) neutrinos emitted by the 
Sun reach the Earth intact and are detected. 

The relative spread in the predictions for the ^Bc neutrino flux <?bc and for 
the signals due to the ''Be neutrinos in the Cl-Ar and Ga-Ge experiments, as it 
follows from Tables 1-3, does not exceed approximately 30%. 

4 The *B Neutrino Problem 

A further inspection of the results collected in Table 1 reveals that the differences 

(and the estimated imccrtaintics) in the predictions of the solar models for the 
total flux of ®B neutrinos are the largest: the value of ^b in the models [29]-[31] is 
by more than a factor of 2 larger than the value obtained in the "low flux" model 
[32]. The ®B neutrinos are born in the Sun in the /J^-decay, ®B ®Be*+e"'" + fc'e, 
of the ®B nucleus which is produced in the reaction 

p+ ^Be^ ®B + 7 (12) 

initiated by ~20 keV protons. Obviously, #b is proportional to the rate of the 
process (12) taking place in the solar plasma environment, which in turn is 
to large extent determined by the cross-section of (12), cri7(Ep). The latter is 
usually represented in the form [23] 

ai7(Ep) = exp(-87r eVv), (13) 

Ep 
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where exp(— Stt e^/v) is the Gamow penetration factor, and Ep and v are the p 
- '^Be cm. kinetic energy and relative velocity. The largely different values of 
the astrophysical factor S17, S17 = Si7(Ep ~ 20 keV), adopted by the authors of 
[29], [30], [31] and of [32] is one of the major sources of the large spread in the 
predictions for #b- 

Because of background problems it is impossible to measure the cross-section 
(Ti7(Ep) directly at the low energies of the incident protons, which are of astro- 
physical interest. The experimental studies of the process (12) were performed 
at energies 110 keV < Ep <; 2000 keV. They are technically rather difficult be- 
cause of the instability of the ''Be serving as a target. The results obtained in 
the indicated higher energy domain are extrapolated to Ep 20 keV using a 
theoretical model describing the data (and the process (12) in the entire en- 
ergy range 20 keV < Ep ^ 2000 keV) and taking into account the possible solar 
plasma screening effects. Obviously, there are at least two major sources of un- 
certainties in the determination of S17 inherent to the indicated approach: the 
uncertainties associated with the data at Ep > 110 keV, and those associated 
with the extrapolation procedure exploited. 

Altogether six experiments have provided data on the the p - ''Be cross- 
section cri7(Ep) so far. The results of the four most accurate of them [44], [45] 
can be grouped in two distinct pairs, [44] and [45], which agree on the energy 
dependence of Si7(Ep), but disagree systematically by ^ (20 25)% ( ^(2-3) 
s.d.) on the absolute values of Si7(Ep). The authors of [29], [31] (and, we suppose, 
of [30]) used in their calculations the value S17 = (22.4 ± 2.1) eV-b derived by 
extrapolation in [46] on the basis of the data from all six experiments. 

A new method of experimental determination of S17 was proposed relatively 
recently in [47]. It is based on the idea of measuring the cross-section of the 
inverse reaction, 7 + ^B ^ ""Be + p, by studying the dissociation of *B into 
p -|- ''Be in the Coulomb field of a heavy nucleus, chosen to be ^"^^Pb. The 
time-reversal symmetry guarantees that the cross-sections of (12) and of the 
inverse reaction should be equal. The extraction of the values of cri7(Ep) (and of 
Si7(Ep)) from the data on the process ^B -|- ^ospb ^ p -|- ^Be-^ 208p^, jg 
straightforward and is associated with certain subtleties (see, e.g., [48]). 

Using the results of the experiment of Motobayashi et al. [47] on the reaction 
83 ^ 208p|3 ^ p ^ '^Be+ 208pb to determine the cross-section £7i7(Ep) in the 
energy interval 500 keV ^ Ep 2000 keV, the results of the most recent of 
the experiments on (12) of Filippone et al. [45] in the interval (110 - 500) keV, 
and a new extrapolation model developed by them, the authors of [32] obtain 
Si7 = 17±2cV-b. " 

The additional difference between the values of <?b predicted in [29] and in 
[32] is due to 



Let us note that in [25] the value of S17 derived in [46] and adopted in [29], [31] 
was also used, but with a larger systematical error, S17 = (22.4 ± 1.3 ± 3.0) eV-b, 
introduced to account for the (systematic) difference between the data on cri7(Ep) 
from the experiments [44] and [45]. 
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i) the use of different (but still compatible within the errors with the measured 
or deduced from the data) values of other relevant nuclear reaction cross- 
sections, as those of the ^He + ^He^ ^He + 2p and ^He + ^He ^ '^Be+7 
reactions (a factor ~ 1.2), and 

ii) the use in [29] and in [32] of different methods to account for the diffusion 
of the heavy elements in the Sun. 

The latter leads, in particular, to a difference in the values of the central tem- 
perature in the Sun in the models [29] and [32]: Tc(BP'95) = 1.584 x 10^ K and 
Tc(DS'96) = 1.561 X lO'^ K. As the neutrino flux is very sensitive to the 
value of Tc, scaling as [40] <?b ~ T^^ , the indicated difference in Tc implies an 
additional difference in the values of predicted in [29] and in [32] (a factor of 
~ 1.4). 

5 The Missing '^Be Neutrinos 

Even if one accepts that there are large uncertainties in the predictions for the 
flux of neutrinos and in all analyses one should rather use for <Pb the value 
implied by the (Super-) Kamiokande data, (2) and (6), another problem arises: 
the predictions of all contemporary solar models for the flux of ^Be neutrinos, 
^Be) considerably larger than the value suggested by the existing solar neu- 
trino data. This was first noticed in [49] and confirmed in several subsequent 
more detailed studies [50]-[53] utilizing a variety of different methods. We shall 
illustrate here the above result using rather simple arguments [49], [53]. 

Let us assume that the spectrum of the ^B neutrino flux coincides with that 
predicted by the solar models, i.e.. with the spectrmn of the i/g emitted in the 
decay ^B ^Be* + e+ -I- z/e (see Fig. 2). This would be the case if the solar 
*B behave conventionally during their journey to the Earth. For the value 
of the total *B neutrino flux one can use the Super-Kamiokande result (6): 
= (2.44 ^0 11) ^ 1*^^ cm^^sec^\ where the statistical and the systematic 
errors were added in quadratures. Knowing <Pb and the the cross-section [23] 
of the Pontecorvo-Davis reaction Vg -1-"^^ CI +'^^ Ar, one can calculate the 

contribution of the ®B neutrinos to the signal in the Davis et al. experiment, 
RB(Ar). One finds: 

R^(Ar) = (2.71 t^oil) SNU. (14) 

By subtracting this value from the rate of Ar production measured in the Davis 
et al. experiment we obtain the sum of the contributions of the ^Be, pep and 
CNO neutrinos to the signal in this experiment: 

j^Be-Fpep-KCNO^^j.) ^ +0.37) gj^y. (15) 

Given the solar model independent relation (9) between ^pcp and ^pp, and 
that #pp is rather tightly constrained by the data on the solar luminosity, we 
can consider as rather reliable (and weakly model dependent) the solar model 
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predictions for Rg^(Ar) = (0.22-0.24) SNU. Taking RP«P(Ar) = 0.22 SNU one 
obtains from (15) 

R^''(Ar) < (-0.37 to.lD SNU, (16) 

R^''(Ar) being the '^Be neutrino contribution to the signal in the Davis et al. 

experiment. At 99.73% C.L. (3 s.d.) this implies R^'=(Ar) < 0.74 SNU, which is 
smaller than the predictions of the solar models with heavy element diffusion, 
RB^(Ar) = (0.89 - 1.24) SNU [28], [29], [30], [31], [32], [35]. As Ri^^(Ar) ~ 
<?Bc, the result obtained suggests that, if the solar neutrinos are assumed to 
behave conventionally on the way to the Earth (i.e., do not undergo oscillations, 
transitions, decays, etc.), the ^Be Ve flux inferred from the solar neutrino data is 
substantially smaller than the flux predicted by the contemporary solar models. 

Similar (though statistically somewhat weaker) conclusions can be reached 
for the contribution of the ^Be neutrinos to the signal in the Ga-Ge detectors, 
R^''(Ge), and correspondingly for <?bc, by taking into account the fact that the 
solar model predictions for the contributions of the pp and pep neutrinos to the 
indicated signal, RPP+P*P(Ge), are tightly constrained by the data on the solar 
luminosity and vary by no more than 3%: Rg^+P*'P(Ge) = (72.7 - 75.0) SNU. 
Subtracting Rg^'''^P(Ge) = 72 SNU from the rate of Ge production observed in 
the SAGE and GALLEX experiments, Eq. (5), one obtains for the contribution 
of **B, ^Bc and CNO neutrinos: RI^+b°+CNO(q^) ^ (3 4 ^ 7) g^y. 

Utilizing the value of measured by the Super-Kamiokande experiment 
and the Ga-Ge reaction cross-section [23], [53] permits to calculate the con- 
tribution of the *^B neutrinos, RB(Ge), to Roxp(Gc): RB(Gc) = (5.9 t^) SNU, 
where the error is dominated by the estimated uncertainty in the value of the 
Ga-Ge reaction cross-section [53]. Subtracting the so derived value of R^(Ge) 
from the value of R^^+^'^+cno (q^) we get: 

R^'(Ge)< (-2.5 t?;^) SNU. (17) 

Consequently, at 99.73% C.L. (3 s.d.) the contribution of ^Be neutrinos to the 
signals in the SAGE and GALLEX experiments docs not exceed 25.1 SNU, while 
the solar models [24]-[32], [36] predict R^^(Ge) > 27 SNU. 

Analogous results have been obtained in [50], [52] using different meth- 
ods. The same conclusion has been reached in [51] as well on the basis of a 
X^— analysis of the solar model description of the data, in which the total pp, 
pep, ^Be, ^B and CNO neutrino fluxes were treated as free parameters subject 
only to the luminosity constraint (11), while the spectra of solar neutrinos were 
assumed to coincide with the predicted ones in the absence of unconventional 
neutrino behaviour (as oscillations in vacuum, etc.). 

Thus, there arc strong indications from the existing solar nciitrino data that 
the flux of ''Be (electron) neutrinos is considerably smaller than the flux pre- 
dicted in all contemporary solar models. Given the results of the GALLEX and 

Actually, the 3 s.d. upper limit on R^''(Ar) is smaller than the predictions of all 
known to the author solar models proposed in the last 10 years (see also the second 
article quoted in [36]). 
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SAGE calibration experiments, we can conclude that both the Davis et al. and 
the Super-Kamiokande (Kamiokande) data, (2) and (6) (Eq. (3)), have to be 
incorrect in order for the above conclusion to be not valid. The discrepancy bet- 
ween the value of suggested by the analyses of the existing solar neutrino 
data and the solar model predictions for <?bo represents the major new aspect 
of the solar neutrino problem. No plausible astrophysical and/or nuclear physics 
explanation of this discrepancy has been proposed so far. 

6 Neutrino Physics Solutions of the Solar Neutrino 
Problem 

We have seen that none of the solar models developed during the last ten years 
provides a satisfactory description of the existing solar neutrino data. The dis- 
crepancy between the data and the solar model predictions is especially large for 
the majority of models with heavy element diffusion, which are compatible with 
the helioseismological data. The solar model predictions for the signals caused 
by the solar neutrinos in the solar neutrino experiments are larger than the mea- 
sured signals. In particular, no solar, atomic or nuclear physics solution to the 
^Be neutrino problem discussed above was found so far. Since the solar neutrino 
detectors are sensitive either only, or predominantly, to the solar Vf^ flux, these 
results indicate that the solar flux is depleted on the way to the Earth. 

Such a depletion can take place naturally if the solar undergo transitions 
into neutrinos of a different type, v^^ and/or v^^ and/or into a sterile neutrino v^^ 
or are converted into antineutrinos v^^ and/or v^-, while they travel to the Earth. 
The depletion of the solar v,, flux might be caused also by instability of the solar 
neutrinos which can decay on their way to the Earth. Thus, several physically 
rather different neutrino physics solutions of the solar neutrino problem are, in 
principle, possible. They all require the existence of "unconventional" intrinsic 
neutrino properties (mass, mixing, magnetic moment) and/or couplings (e.g., 
flavour changing neutral current (FCNC) interactions). More speciflcally, these 
solutions include: 

i) oscillations in vacuum [9], [16], [2] of the solar into different weak eigen- 
state neutrinos {y^ and/or Vr^ and/or sterile neutrinos, on the way from 
the surface of the Sun to the Earth [54]-[58], 

ii) matter-enhanced transitions [38], [39] i/g '^h(t)j and/or i/g ^ i/^, while the 
solar neutrinos propagate from the central part to the surface of the Sun 
[59]-[61], 

iii) solar Ue resonant spin or spin-flavour conversion (RSFC) in the magnetic 
field of the Sun [62] , and 

iv) matter-enhanced transitions, for instance Ue — * ^^t, in the Sun, induced by 
flavour changing neutral current (FCNC) interactions of the solar Ve with 
the particles forming the solar matter [63]-[65] (these transitions can take 
place even in the case of absence of lepton mixing in vacuum and massless 
neutrinos [63]). 
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All these possibilities have been and continue to be extensively studied (see the 
quoted articles). There have not been recent studies of the solar neutrino decay 
hypothesis [66] which, however, was disfavored [67] by the earlier solar neutrino 
data. 

In what follows we shall discuss the vacuum oscillation and the matter- 
enhanced transition solutions of the solar neutrino problem. The status of these 
solutions has been reviewed recently, e.g., in [68], [69]. 

6.1 Oscillations in Vacuum 

Neutrino oscillations in vacuum [9] , have been discussed in connection with the 
solar neutrino experiments [16] and as a possible solution of the solar neutrino 
problem [2], [1], [54]-[58], [69] (and the literature quoted therein) for about 31 
years. In the simplest version of this scenario it is assumed that the state vector 
of the electron neutrino, jj/g), produced in vacuum with momentum p in some 
weak interaction process, is a coherent superposition of the state vectors jfi) of 
two neutrinos Vi, i=l,2, having the same momentum p and definite but difi'erent 
masses in vacuum, mj, mi ^ m2, while the linear combination of and [^2), 
which is orthogonal to ji'e), represents the state vector \ux) of another weak- 
eigenstate neutrino, l^x) = Wh{t)) ot \i/s)' 

= cos^ + [^'2) sin^ , (18a) 

= — sin^ + [1/2) cos^ , (18&) 

where Q is the neutrino (lepton) mixing angle in vacuum. We shall assume for 
concreteness in what follows that is an active neutrino, say f^, = 

Obviously, the states |i'i,2) are eigenstates of the Hamiltonian of the neutrino 
system in vacuum, Hq: 

Ho\vi) =Ei\ui), Ei = ^p^ + ml i = l,2. (19) 

If Vf, is produced at time f = in the Sun in the state given by (18a), after 
a time t the latter will evolve into the state 

\iy^{t)) = e-'^i* ]z/i) COS 61 + e"'^^* jz/a) sin 61 , (20) 

where we have ignored the overall space coordinate dependent factor exp(ipr) 

in the right-hand side of (20) and have assumed that the solar matter docs not 
affect the evolution of the neutrino system. (The possible effects of matter on 
the evolution of the neutrino state will be considered in the next Section.) Using 
(18a) and (18b) to express the vectors ji^i) and [1/2) in terms of the vectors Ivg) 
and jfju) we can rewrite (20) in the form: 



where 



Aee(i) = e-*-^i* cos^ e + e"*-^^* sin^ 9 



(22) 
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and 

= I sin 26 (e-*^=* - e"*^^*) (23) 

arc the probability amplitudes to find respectively neutrino i/e and neutrino 
at time t of the evolution of the neutrino system if neutrino u^. has been produced 
at time t = 0. Thus, if neutrinos z^i and 1/2 are not mass-degenerate, mi 7^ m2, 
and if nontrivial neutrino mixing exists in vacmim, 6 7^ mr/2, n = 0, 1, 2, we 
have |Ajue(t)|^ ^ and transitions in flight between the states Ivg) and li^^) (i.e., 
between the neutrinos and f^u) are possible. 

Assuming that neutrinos i>i and 1^2 are stable and relativistic, it is not difficult 
to derive from (22) and (23) the probabilities that a solar i/g with energy E = 
\p\ = p will not change into on its way to the Earth, Pyoi'^e ~* ^e:t)j and 
will transform into while traveling to the Earth, Pyoi'^e — * i^n',t): 

Pvo{ye ye;t) = \A^e{t)\^ = 1 -^s.m^26 

Pvoi^e ^ t) = |A^e(i)|^ = \ sin^ 2^ 
where Am^ = m\ — m\, 

Ly = A-kE/ Am- 
is the oscillation length in vacuum, 

R{ty) = i?o [1 - ecos(27rtj,/T)] , (27) 

is the Sun-Earth distance at time ty of the year (T = 365 days), i?o = 1.4966 x 10^ 
km and e = 0.0167 being the mean Sun-Earth distance and the elUpticity of the 
Earth orbit around the Sun. In deriving (24) and (25) we have used the equalities 

E2-Ei^p + Am^/{2p) and t ^ R{ty) 

valid for relativistic neutrinos 2^1,2. The quantities Am? and sin^ 26 are typically 
considered and treated as free parameters to be determined by the analysis of 
the solar neutrino data. 

It should be clear from the above discussion that the neutrino oscillations, 
if they exist, would be a purely quantum mechanical phenomenon. The require- 
ments of coherence between the states \v\) and \u2) in the superposition (18a) 
representing the p,. at the production point, and that the coherence be main- 
tained during the evolution of the neutrino system up to the moment of neutrino 
detection, are crucial for the neutrino oscillations to occur. The subtleties and 
the implications of the coherence condition for neutrino oscillations continue to 
be discussed (see, e.g., [2], [70], [71] and the articles quoted therein). 

As it follows from (25), the Ue transition probability Pvo{ve i^fi'tt), 

depends on two factors: on (1 — coa2TrR{ty)/ Ly), which exhibits oscillatory de- 
pendence on the distance traveled by the neutrinos and on the neutrino energy 
(hence the name "neutrino oscillations"), and on sin^ 2^ which determines the 



1 — cos 27r 



R{ty) 



(24) 



l_cos27r:^ 



(25) 



(26) 
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amplitude of the oscillations. In order for the Vfj, oscillation probability 

to be large, Pvo{i'e ^ii]t) — !> two conditions have to be fulfilled: the neu- 
trino mixing in vacuum must be large, sin^ 2^ = 1, and the oscillation length 
in vacuum L^. has to be of the order of or smaller than the distance traveled 
by the neutrinos, R: <> 2nR. If the second condition is not satisfied, i.e., if 
Ly » 27ri?, the oscillations do not have enough time to develop on the way to the 
neutrino detector as the source-detector distance R (in our case the Sun-Earth 
distance) is too short, and one has Pvoit^e i^n'^t) — 0. 

Let us note that, in general, a given experiment searching for neutrino oscil- 
lations, is specified, in particular, by the average energy of the neutrinos being 
studied, E, and by the distance traveled by the neutrinos to the neutrino detec- 
tor. The requirement Ly ^ 2-kR determines the minimal value of the parameter 
Am^ to which the experiment is sensitive (figure of merit of the experiment): 
min(Zim^) ~ 2E/R. Because of the interference nature of the neutrino oscilla- 
tions, the neutrino oscillation experiments can probe, in general, rather small 
values of Am? (see, e.g., [1], [2]). In addition, due to the large distance be- 
tween the Sun and Earth and the relatively low energies of the solar neutrinos, 
~ 1 MeV, the experiments with solar neutrinos have a remarkable sensitivity 
to the parameter Am?, namely, they can probe values as small as 10~^^ eV^: 
Am"^ ^ 10-" eV2. 

To summarize the above discussion, if (18a) is realized and Am^ ^ 10"^^ eV^ 
the solar can take part in vacuum oscillations on the way to the Earth. In 
this case the flavour content of the electron neutrino state vector will change 
periodically between the Sun and the Earth due to the different time evolution 
of the vector's massive neutrino components. The amplitude of these oscillations 
is determined by the value of sin^ 29. If sin^ 29 is sufficiently large, the neutrinos 
that are being detected in the solar neutrino detectors on Earth will be in states 
representing, in general, certain superpositions of the states of^ Ve and Vf^. As 
the muon (and tau and sterile) neutrinos interact much weaker with matter than 
electron neutrinos, the measured signals in the solar neutrino detectors should 
be depleted with respect to the expected ones. This would explain the solar 
neutrino problem. 

Detailed analyses of the solar neutrino data in terms of the hypothesis of 
two-neutrino vacuum oscillations of solar neutrinos have been performed in the 
period after 1991, e.g., in [54], [55], [57], [68], [69]. It was found that the two- 
neutrino oscillations involving the Ug and an active neutrino, Ve ^'/^(t); provide 
a good quality description (x^-fit) of the solar neutrino data for values of the 
two vacuum oscillation parameters belonging approximately to the region (see, 
e.g., [68]): 

5.0 X lO-^eV^ Am^ <> 10"^^, (28a) 

0.65^ sin^ 261 < 1.0. (286) 

" Obviously, if i>e mixes with Vfj, and /or Vr and/or Vs, these states will be superpositions 
of the states of f^, and/or Vr and/or i/g. 
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Energy, MeV 



Fig. 5. The vacuum oscillation probability Pvoife — > fe'^t), Eq. (24), for the mean 
distance between the Sun and the Earth, t = (upper frame), and the probability 
Pvoi^'e i^e'jt) averaged over a period of 1 year (lower frame), as functions of the 
neutrino energy E for Arn^ = 10"^° eV^ and sin^ 261 = 0.8 (from [57]). For further 
details see the text. 

At the same time, as it was shown in [55], [57], the oscillations into sterile neu- 
trino Us, Ve t/g, give a poor fit of the solar neutrino data and are thus strongly 
disfavored by the data as a possible solution of the solar neutrino problem. 

The probability of solar Ue survival, Pyoi^e i^e',Ro), in which t = R{ty) is 
replaced with the average Sun- Earth distance Rq, and the probability Pyoi'^e 
V(.: t) averaged over the period of one ycar,-'^ arc shown for Am^ = 6.3 x lO^^^cV^ 
and sin^ 29 = 0.8 as functions of the solar neutrino energy E in the upper and 
lower frames of Fig. 5, respectively (taken from [57]). 

Although in the analyses [57] , [68] leading to the above results the predictions 
of the solar model [29] with heavy element diffusion for the fluxes of the pp, pep, 
^Be, and CNO neutrinos were used, it was also verified [57], [58], [68] that 
the results so obtained (i.e., the existence of the vacuum Uf, ^ i^h{t) oscillation 

^ The one-year averaged probability has to be used in the analyses of data taken over 
periods of k years, k = 1,2,3,..., as are the data (1) - (4) provided by the Cl-Ar, 
Ga-Ge and Kamiokande experiments. 
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^^^1 4 experiments + Super-Kamiokande 

(exp/SSM) = 0.400±0.024 

SSM: Bahcall and Pisonneault 1995 
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Fig. 6. Regions of values of the parameters Am^ and sin^ 29 (shown in black) for which 
the solar neutrino data can be described (at 95% C.L.) in terms of vacuum ^ i'/jCt) 
oscillations of the solar Ve (from [72]). For further details see the text. 



solution) are stable with respect to variations of the values of the total and ^Be 
neutrino fluxes within wide intervals in which the predictions of all contemporary 
solar models lie.^ 

The results of a recent analysis of the solar neutrino data in terms of the 
hypothesis of two-neutrino Ve ^ ^fj-ir) oscillations of the solar neutrinos [72], are 
shown graphically in Fig. 6. The analysis was based on the predictions of the 
solar model of Bahcall and Pinsonneault [29] with heavy element diffusion. The 
regions in the Am^ — sin^ 29 plane colored in black correspond to values of the 
two parameters for which one obtains (at 95% C.L.) a description of the data, 



" Actually, in the analysis performed in [57] the *B neutrino flux <Pb was treated as a 
free parameter, while the ^Be neutrino flux ^bc was assumed to take values in the 
interval (0.7 - 1.3)<5i^^^ where <1>W'^^ is the flux in the model [29] (see Table 1). 
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in other words, a solution of the solar neutrino problem.'' 

Let us note that for the values of Am"^ from the interval (28a), the oscillation 
length in vacuum for the solar neutrinos with energy E 1 MeV is of the order 
of the Sun-Earth distance: ~ (2.5 — 5.0) x lO'^ km. At the same time Ly is 
much bigger than the solar radius: Ly ^ Rq. 

6.2 Matter-Enhanced Transitions 

Let us consider next that the possible effects of the solar matter on the oscilla- 
tions of solar neutrinos assuming that (18a) and (18b) hold true and supposing 
first that \iyx) = W^i)- 

The presence of matter can drastically change the pattern of neutrino oscilla- 
tions: neutrinos can interact with the particles forming the matter. Accordingly, 
the Hamiltonian of the neutrino system in matter differs from the Hamiltonian 
of the neutrino system in vacuum Hq, 

Hm = Ho + Hint , (29) 

where Hint describes the interaction of the flavour neutrinos with the particles 
of matter. When, e.g., electron neutrinos propagate in matter, they can scat- 
ter (due to the Hint) on the particles present in matter: on the electrons (e~), 
protons {p) and neutrons (n). The incoherent elastic and the quasi-elastic scat- 
tering, in which the states of the initial particles change in the scattering process 
(destroying the coherence between the neutrino states), are not of interest for 
our discussion for one simple reason - they have a negligible effect on the solar 
neutrino propagation in the Sun * : oven in the center of the Sun, where the den- 
sity of matter is relatively high 150 g/cuv^), an I'e with energy of 1 MeV has a 
mean free path with respect to the indicated scattering processes, which exceeds 
10^'^ km (recall that the solar radius is much smaller: Rq = 6.96 x 10^ km). The 
oscillating Uf. and can scatter also elastically in the forward direction on the 
e~, p and n, with the momenta and the spin states of the particles participating 
in the elastic scattering reaction remaining unchanged. In such a process the 
coherence of the neutrino states is being preserved and the oscillations between 
the flavour neutrinos can continue in spite of, and actually, in parallel to, the 
scattering. 

The Uf. and coherent elastic scattering in the forward direction on the 
particles of matter generates nontrivial indices of refraction of the Ve and in 
matter [38]: K(fe) 7^ 1, 1- Most importantly, the index of refraction of 

the Ve thus generated does not coincide with the index of refraction of the f^: 

In this analysis the experimental results (1), (2), (4) and (6) and the data from 
the GALLEX experiment obtained in 51 runs of measurements, Rgallex(Go) = 
(69.7 ±6.7l ID SNU, were used. The solution regions of values of Am^ and sin^ 29 
do not change significantly if one uses the most recent GALLEX data, (4). 
' These processes are important, however, for the supernova neutrinos (see, e.g., [5], 
[6]). 
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K(i^e) 7^ ^(f^). The difference between the two indices of refraction is determined 
essentially by the difference of the real parts of the forward Ve — e~ and — e~ 
elastic scattering amplitudes [38], Re [Fj,^_e-(0)] — Re [F^^_e-(0)]: due to the 
flavour symmetry of the neutrino - quark (neutrino - nucleon) neutral current 
interaction, the forward Ug — p,n and — p,n elastic scattering amplitudes 
are equal and therefore do not contribute to the difference of interest.-' The 
real parts of the amplitudes Fj^^__e-(0) and F^^_q-{0) can be calculated in the 
standard theory. One finds the following result [38], [73], [74] (see also [75]) for 
the difference of the indices of refraction of Ve and v^: 

2tt 1 
K(z^e) - k(z^^) = ^ (Re [F,,_e- (0)] - Re [F,^_e- (0)]) = - -V2GFNe , (30) 

where Gp is the Fermi constant and Ng is the electron number density in 

matter. Let us note that the forward scattering amplitudes for the antineu- 
trinos Fp^_e-(0) and Fp^_g-(0) coincide in absolute value with the amplitudes 
F^^_e-(0) and Fy^_e-(0) but have opposite sign and therefore one has 

«(i>e) - K(i^/*) = + ^y^Gp Ne . (31) 

Knowing the expression for the difference of the indices of refraction of Ve 
and in matter, it is not difficult to write the system of evolution equations 
which describes the Ve ^ oscillations in matter [38], [73], [74]: 

d ( A,{t, to)\_[ -e{t) e' \f Ae{t, h) \ . . 



where Ae{t, to) {Af^{t, to)) is the amplitude of the probability to find neutrino Ve 
{u/j,) at time t of the evolution of the neutrino system if at time to the neutrino 
Ve or (or a state representing a linear combination of the states describing 
the two neutrinos) has been produced, t > to. Furthermore e{t) and e' are real 
functions of the neutrino cnc;rgy E = p. of Am'^ , of the mixing angle in vacuum 
9 and of the electron number density in the point of the neutrino trajectory in 
matter reached at time t, Ne{t), 

e{t) = -[^cos2e-V2GFNe{t)], e' = ^sm20. (33) 



We standardly assume that the weak interaction of tlic flavour neutrinos Ue, i^n and 
Vt and antineutrinos Ve, i>fj. and Vt is described by the standard (Glasiiow-Salam- 
Weinberg) theory of electroweak interaction and that the generation of nonzero neu- 
trino masses and lepton mixing leading to (18a) and (18b) does not produce new 
couplings which can change substantially the neutrino weak interaction, as required 
by the existing experimental limits on such new couphngs (for an alternative possi- 
bility see, e.g., [63]). Let us add that the imaginary parts of the forward scattering 
amplitudes (responsible, in particular, for decoherence effects) are proportional to the 
corresponding total scattering cross-sections and in the case of interest are negligible 
in comparison with the real parts. 
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The term \/2GFNf,{t) in the parameter e{t) accounts for the effects of matter 
on the neutrino oscillations. 

Let us note that the system of evolution equations describing the oscillations 
of antineutrinos f'/i in matter has exactly the same form except for the 

matter term in e{t) which, in accordance with (30) and (31), changes sign. 

Due to the presence of the interaction term Hint in the Hamiltonian of the 
neutrino system in matter Hm , the eigenstates of the Hamiltonian of the neutrino 
system in vacuum, and 11^2), are not eigenstates of Hm- As a result of the 
coherent scattering of and off the particles forming the matter transitions 
between the states \ui) and \u2) become possible in matter: 

{l^2\ Hint Wl) 7^ 0. (34) 

Consider first the case of Ue <-> oscillations taking place in matter with electron 
number density which does not change along the neutrino trajectory: Ne{t) = 
Ne = const. It proves convenient to find the states |!^™2)j which diagonalize 
the evolution matrix in the right-hand side of the system (32), or equivalently, 
the Hamiltonian of the neutrino system in matter. The relations between the 
matter-eigenstates |j^™2) ^-nd the flavour-eigenstates l^cij,) have the same form 
as the relations (18a) and (18b) between the vacuum mass-eigenstates |i^i,2) and 

We) = WT)cos0m + WT)sin9m, (35a) 
W^^) = -WD sin^m + WT) cos^™ . (356) 
Here 6m is the neutrino mixing angle in matter [38], 
. „n e' tan 2^ 

sin2g^= , ^ , (36) 

^(l-^)2+tan^2^' 

cos2^„ = -.=L= = " (37) 



^(l-^)2+tan2 2^' 



where the quantity 



2EV2Gf 

is called "resonance density" [73]. The mattcr-cigcnstatcs |z^™2) (which arc also 
called "adiabatic") are eigenstates of the evolution matrix (Hamiltonian) in (32), 
corresponding to the two eigenvalues, i^i2, whose difference is given by 



- ET = 2V^T7^ = ^ ^(1 _ _^)2 eos2 26 + sin' 26. (39) 

It should be almost obvious from (25) after comparing (18a), (18b) with 
(35a), (35b) that the probability to find neutrino at time t if neutrino has 
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been produced at time to = and it traversed a distance {t — to) = t = Rm in 
matter with constant electron number density A^e, has the form [39]: 



1 o 

1 — COSZTT — 



(40) 



where 

27r^ = {E^ - ET)t, (41) 



and is the oscillation length in matter: 



(1- j^)2cos2 26' + sin2 261 



(42) 



Evidently, the amplitude of the <-> oscillations in matter is equal to 
sin^ 29m- It follows from (36) that, most remarkably, the dependence of sin^ 29jn 
on Ne has a resonance character [39]. Indeed, if in the case of interest the con- 
dition 

Am^ cos^ 29 >0 (43) 

is fulfilled, for any finite value of sin^ 26 there exists a value of the electron 
number density equal to N^^^, such that when 

Ne = Nl""^ (44) 

we have 

sin^ 26m = l- (45) 

Note that if A^e = A^J^* , wc get sin^ 29„i ~ 1 even if the mixing angle in vacuum 
is small, i.e., if sin^ 20 <C 1. This implies that the presence of matter can lead to 
a strong enhancement of the oscillation probability Pm{i^e — * t^ti'^t) even when 
the ^ i^p, oscillations in vacuum are strongly suppressed due to a small value 
of sin^ 29 (hence, the name "matter-enhanced neutrino oscillations" ) . 
The oscillation length at resonance is given by 



sin 20 



while the width in TV,, of the resonance (i.e., the "distance" in Ng between the 
points at which sin^ 29m = 1/2) has the form 

^j^rea ^ tan 26 . (47) 

Thus, if the mixing angle in vacuum is small the resonance is narrow, AN^'^'^ <^ 
N^'^", and the oscillation length in matter at resonance is relatively large, L^" ^ 
Ly. As it follows from (39), the energy difference — £™ has a minimum at 
the resonance: 

{E^ - EY'Y''' = min {E^ - £]") = sin 20. (48) 
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It is instructive to consider two limiting case. If Ne <C N^^^ , as it follows from 
(36) and (42), 9m — (sin20m — sin2^^), = Lv and the neutrinos oscillate 
practically as in vacuum. In the opposite limit, A^e ^ A^J'''*, iV™*tan^20, one 
finds from (36) and (37) that 6^ = 1^/2 (sin2^?„J = 0, cos2^?,„ = —1) and the 
presence of matter suppresses the i/g ^ '^/j. oscillations (see (40)). In this case we 
get from (35a) and (35b): 

We) = WT), (49a) 

i.e., if the electron number density exceeds considerably the resonance density, 
Ve practically coincides with the heavier of the two matter-eigenstate neutrinos 
f™, while the coincides with the lighter one z/J". 

The analogs of (36), (37), (39), (40) and (42) for oscillations of antineutrinos, 
f'e f'/x' matter with constant A^e can formally be obtained by replacing A^e 
with {—Ne) in the indicated equations. If condition (43) is fulfilled, we have 
N^''^ > and the term (1 + Ne/Nl^^) which appears, e.g., in the expression for 
the mixing angle in matter Om in the case of <-> C*^ oscillations, can never be 
zero. Thus, a resonance enhancement of the T>e f'/i oscillations cannot take 
place. The matter, actually, can only suppress the oscillations. 

It should be clear from this discussion that depending on the sign of the 
product Am^ cos 20, the presence of matter can lead to resonance enhancement 
either of the or of the Dg v^^ oscillations, but not of the both types 

of oscillations. This is a consequence of the fact [76] that the matter in the Sun 
or in the Earth we are interested in, is not charge-symmetric (it contains e~ , p 
and n, but does not contain their antiparticles) and therefore the oscillations in 
matter are neither CP- nor CPT- invariant.*^ In what follows we shall assume 
that A-m? > and cos 20 > 0, so that (43) is satisfied and therefore only the 
t'e i^n oscillations can be enhanced by the matter effects. 

Since the neutral current weak interaction of neutrinos in the standard theory 
is flavour symmetric, the formulae and results we have obtained above and shall 
obtain in what follows arc valid for the case of i^g — mixing ((18a) and (18b)) 
and Ve i^T oscillations in matter as well. In what concerns the possibility 
of mixing and oscillations between the and a sterile neutrino Vg, i^e t^s, 
the relevant formulae can be obtained from the formulae derived for the case of 
fe '^nir) oscillations by [76] replacing with (A'^e — l/2A^n)! where A^^^ is the 
number density of neutrons in matter. 

The formalism wc have developed above can be directly applied, for instance, 
to the study of the matter effects in the '^ij.(t) {'^ij.(t) ^ ^e) oscillations of 

the flavour neutrinos which traverse the Earth mantle (but do not traverse the 
Earth core). The electron number density changes little around the mean value 

* As it is not difEcult to convince oneself, the matter effects in the Ve ^ {ve ^ v^) 
oscillations will be invariant with respect to the operation of time reversal if the A^e 
distribution along the neutrino path is symmetric with respect to this operation. 
The latter condition is fulfilled for the Ne distribution along a path of a neutrino 
crossing the Earth [77]. 
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of -/Ve = 2.3 cm~^ Na, Na being the Avogadro number, along the trajectories 
of neutrinos which cross a substantial part of the Earth mantle and the Nf. = 
const, approximation is rather accurate. If, for example. Am? = 10~^ eV^, E = 
1 GcV and sin^ 261 ^ 0.5, wc have: N^"-' ^ 4.6 cmr^ Na. sin^ 26^ = 0.8 and the 
oscillation length in matter, L„ = 3 x 10^ km, is of the order of the depth of the 
Earth mantle,' so that one can have ^wRm ^ L^- It is not difficult to obtain 
an expression for the ly^ ^ '^/j(t) oscillation probability in the case when the 
neutrinos traverse both the Earth mantle and the core assuming A^e is constant, 
but has different values in the two Earth density structures. 

It is not clear, however, what the above interesting results have to do with the 
problem of main interest for us, namely, accounting for the effects of solar matter 
in the oscillations of solar neutrinos while they propagate from the central part to 
the surface of the Sun. The electron number density (the matter density) changes 
considerably along the neutrino path in the Sun: it decreases monotonically from 
the value of ~ 100 cm~^ Na (^ 150 g/cm.^^) in the center of the Sun to at the 
surface of the Sun. Actually, according to the contemporary solar models (see, 
e.g., [23], [29]), Nf. decreases approximately exponentially in the radial direction 
towards the surface of the Sun: 



where {t — to) = cZ is the distance traveled by the neutrino in the Sun, Ae(io) is 
the electron number density in the point of neutrino production in the Sun, ro 
is the scale- height of the change of Ne{t) and one has [23] rg ~ O.li?0. 

Obviously, if Ng changes with t (or equivalently with the distance) along the 
neutrino trajectory, the matter-eigenstates, their energies, the mixing angle and 
the oscillation length in matter, become, through their dependence on Ag, also 



functions of t: [z^i^z) = KW)), = ^^^(O, = 0^{t) and L„ = L^{t). 



It is not difficult to understand qualitatively the possible behaviour of the 

neutrino system when solar neutrinos propagate from the center to the surface 
of the Sun if one realizes that one is dealing effectively with a two-level system 
whose Hamiltonian depends on time and admits "jumps" from one level to the 
other (sec (32)). Let us assume first for simplicity that the electron number 
density in the point of a solar v^. production in the Sun is much bigger than the 
resonance density, Nf.{to) A^^*, and that the mixing angle in vacuum is small, 
sm9 ^ 1. Actually, this is one of the cases relevant to the solar neutrinos. In this 
case we have Omito) = 7r/2 and the state of the electron neutrino in the initial 
moment of the evolution of the system practically coincides with the heavier of 
the two matter-eigenstates: 




(50) 




i^^ito)). 



(51) 



The Earth radius is 6371 km; the Earth core, whose density (Ae) is larger approxi- 
mately by a factor of 2.5 than the density (Ag) in the mantle, has a radius of 3486 
km, so the Earth mantle depth is 2885 km. 
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Thus, at to the neutrino system is in a state corresponding to the "level" with 
energy E^{to). When neutrinos propagate to the surface of the Sun they cross 
a layer of matter in which Ne = Nl"^: in this layer the difference between the 
energies of the two "levels" [El]^ {t) — [t]) has a minimal value on the neutrino 
trajectory ( (39) and (40)). Correspondingly, the evolution of the neutrino system 
can proceed basically in two ways. First, the system can stay on the "level" with 
energy E™{t), i.e., can continue to bo in the state |i^™(t)) up to the final moment 
ts , when the neutrino reaches the surface of the Sun. At the surface of the Sun 
N,{ts) = and therefore Omits) = 0, |j^i™2(is)) = 1^^1,2) and EY:2{ts) = £^1,2- 
Thus, in this case the state describing the neutrino system at to will evolve 
continuously into the state 11^2) at the surface of the Sun. Using (18a) and (18b), 
it is trivial to obtain now the probabilities to find respectively neutrino Vg and 
neutrino at the surface of the; Sun (given the fact that has been produced 
in the initial point of the neutrino trajectory): 

P{ue u^]tsM) = \A^{ts,tQ)? = \WeW2)f = sm^O, (52a) 

P{ue ^ u^;ts,to) = \A^{ts,to)\'' ^ \W^,W2)\'' = cosH. (526) 

It is clear that under the assumptions made (i.e., sin^ ^ ^ 1), a practically total 

Ve — conversion is possible in the case under study. This type of evolution of 
the neutrino system as well as the v,, transitions taking place during the 

evolution, are called [39] "adiabatic". They are characterized by the fact that 
the probability of the "jump" from the upper "level" (having energy E'^{t)) to 
the lower "level" (with energy E^{t)), P' , or equivalently the probability of the 
i^r(to) Kits) transition, P' = P'(i/^(to) ^ t^^its)), on the whole neutrino 
trajectory is negligible: 

P' = P'iv^ito) ^ vTits)) =0 ■■ adiabatic transitions. (53) 

The second possibility is realized if in the resonance region, where the two 
"levels" approach each other most (the difference between the energies of the 
two "levels" iE^it) — E^it)) has a minimal value), the system "jumps" from 
the upper "level" to the lower "level" and after that continues to be in the state 
|i/™(t)) until the neutrino reaches the surface of the Sun. Evidently, now we have 
P' = P'iv^ito) i/™(ts)) = 1. In this case the neutrino system ends up in the 
state |z^i"(ts)) = at the surface of the Sun and the probabilities to find the 
neutrinos and at the surface of the Sun are given by 

P(z/e ^ iye;ts,to) = |Ae(ts,to)|^ = \{VeWl)\'^ = COs'^ 0, (54a) 

Piue ^ u^;ts,to) = |A^(t„to)r = K^^mI^i)!" = sin^^- (546) 

Obviously, if sin^ ^ ^ 1, practically no transitions of the solar i/g into i/^ will 
occur. The considered regime of evolution of the neutrino system and the cor- 
responding Uf. transitions are usually referred to as "extremely nonadia- 
batic". 
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Clearly, the value of the "jump" probability P' plays a crucial role in the the 
— * transitions: it fixes the type of the transition and determines to large 
extent transition probability. We have considered above two limiting 

cases: P' = and P' = 1. Obviously, there exists a whole spectrum of possibilities 
since P' can have any value from to 1. In general, the transitions are called 
"nonadiabatic" if P' is non- negligible (see further). 

Numerical studies have shown [39] that solar neutrinos can undergo both 
adiabatic and nonadiabatic i/g transitions in the Sun and the matter 

effects can be substantial in the solar neutrino oscillations for a remarkably wide 
range of values of the two parameters Aw? and sin^ 20, namely for 

10-^eV^ ^ Z\m^ ^ 10-*eV^ (55a) 

10"* <, sin^ 261 < 1.0 . (556) 

It would be preferable to make more quantitative the preceding analysis. We 
will obtain first the adiabaticity condition [39], [78]. 

Using the (35a) and (35b) we can express the probability amplitudes Ae{t, to) 
and An{t, to) in terms of the probability amplitudes Ai{t, to) and A2{t, to) to find 
the neutrino system in the states |z^i"(t)) and [^'^(i))), respectively, at time t: 

Ae{t,to) = Ai{t,to)cosernit)+A2{t,to)smernit), (56a) 

Ai,{t,to) = -Ai{t,to)smem{t)+A2it,to)cosem{t). (566) 

Substituting (56a) and (56b) in (32) we obtain the system of evolution equations 
for the probability amplitudes Ai{t,to) and ^2(^,^0): 



I (A,{t,to)\ ^ (Er{t) -iO^m (A,{t,to)\ , . 

dt \A2{t,to)J \iem{t) E^{t) ) \A2{t,to)J • ^ ' 

Here ()m{t) = -^0m{t)- It follows from the preceding discussion that the solar neu- 
trino transitions in the Sun will be adiabatic (nonadiabatic) if the nondiagonal 
term in the evolution matrix in the right-hand side of (57), which is responsible 
for the i'2^{to) — > z^r'(is) transitions, is sufficiently small (is non- negligible) . The 
corresponding conditions can be written as 

4n(f) 1, adiabatic transitions, (58a) 

4n(f) ^ 1, nonadiabatic transitions, (586) 
where the adiabaticity function 4n(t) is given by 

4n(0 = ^^^^'^ = V2Gf^^S^ tan^ 29 (l + tan-^20^{t))^ . (59) 

^ 2\em{t)\ mt)\ ^ "^^^ 

In (59) Ne{t) = ffNeit) and we have used (36), (37) and (39) to derive it. 
Expression (59) for 4n{t) implies that the solar neutrino transitions in the Sun 
will be adiabatic if the electron number density changes sufiiciently slowly along 
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the neutrino trajectory; if the change of Ne{t) is relatively fast, the transitions 
would be nonadiabatic. 

In order for the solar neutrino transitions to be, e.g., adiabatic, condition 
(58a) has to be fulfilled in any point of neutrino trajectory in the Sun. However, 
it is not difficult to convince oneself using (36), (37), (50) and (59) that if the 
solar neutrinos cross a layer with resonance density N^^^ on their way to the 
surface of the Sun, condition (58a) will hold if it holds at the resonance point, 
i.e., for the parameter 



rj\Tres\2 

4no = 4n(t = Ues) = V2Gf . " ^ tan^ 261 

mt = tres)\ 

Am^ sin^ 26» _ Ar''^'' 
~ IE cos 26* ~ ^ ' 



where t^es is the time at which the resonance layer is crossed by the neutrinos, 
to < tres < ts, Ar"""" = 2{Nl'"' /\Ne{t = tres)\) tan 20 = 2ro tan 26* is the spatial 
width of the resonance and we have used (38) and (46). Thus, the value of the 
adiabaticity parameter An^ determines the type of the solar neutrino transitions. 
It follows from (60), in particular, that the transitions will be adiabatic if the 
width of the resonance is bigger than the oscillation length at resonance. 

Actually, the system of evolution equations (32) can be solved exactly in 
the case when Ne changes exponentially, (50), along the neutrino path in the 
Sun [79], [80]. On the basis of the exact solution, which is expressed in terms 
of confluent hypergeometric functions [81], it was possible to derive a complete, 
simple and very accurate analytic description of the matter-enhanced transitions 
of solar neutrinos in the Sun [79], [82]-[85] (for a review see [86]). The probability 
that a Ve having momentum p (or energy E = p) and produced at time to in the 
central part of the Sun will not transform into t'^(r) on its way to the surface of 
the Sun (reached at time tg) is given by 

^©(fc'e ^ fc'e; ^s) io) = -P©(^'e ^ i^e; is) io) + Oscillating terms. (61) 

Here 

pQiUe ^ l^e;ts,to) = -P© = ^ + - P') COS 2^„ (to) COS 2^ (62) 

is the average probability, where 



exp 



-27rro 



2E 



exp 



-27rro 



Am-' 



2E 



1 — exp [— 27rro 



2E 



exp [— 27rno(l — tan^ 6)\ — exp [— 27rno(tan 6 — t&v? 6)\ 
1 - exp [-27rno(tan~^ 6 - tan^ 6)\ 



(63) 
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is [79], [86] the "jump" probability for exponentially varying electron number 
density"* A^e, and 9m{to) is the neutrino mixing angle in matter in the point of 
Ve production. 

Wc will not give the explicit analytic expressions for the oscillating terms 
in the probability P(v;{ve z^ei^sj^o)) although they have been derived in the 
exponential density approximation for the N^, as well [84] (see also [89] ) . These 
terms were shown [85] to be, in general, strongly suppressed by the various 
averagings one has to perform when analyzing the solar neutrino data in terms 
of the hypothesis that solar neutrinos undergo matter-enhanced transitions in 
the Sun. More specifically, it was found [85] that the oscillating terms in PQiiu,, 
Ue'itsj to) can be important only for the monochromatic ^Be- and pep-neutrinos 
and only for values of Am^ ^ 10~* eV^. As we shall see, the current solar 
neutrino data suggest that Z\m^ 10^''' eV^. 

It should be emphasized that for Am^ <; 10"'' eV^ the averaging over the 
region of solar neutrino production in the Sun and the integration over the 
neutrino energy renders negligible all interference terms which appear in the 
probability of survival due to the ^ ^h{t) oscillations in vacuum taking 
place on the way of the neutrinos from the surface of the Sun to the surface 
of the Earth. Thus, the probability that will remain while it travels from 
the central part of the Sun to the surface of the Earth is effectively equal to 
the probability of survival of the Vf, while it propagates from the central part 
of the Sun to the surface of the Sun and is given by the average probability 
pQiye '^e',ts,to) (determined by (62) and (63)). 

The probability pQive — * ^e'^ts^to) has several interesting properties. If the 
solar Ve transitions are adiabatic (i.e., P' ^ 0) and cos20„i(to) — —1 (i.e., 
N(.{to)/N^''^ ^ 1, tan 20, solar neutrinos are born "above" and "far" (in A^e) 
from the resonance region), one has 

Pive^Ve;ts,to)=sm'^e, (64) 

which is compatible with the qualitative result (52a) derived earlier. The solar i/^ 
undergo extreme nonadiabatic transitions in the Sun (4no 4C 1) if, e.g., E/Aw? 

is "large" (see (60)). In this case again cos26'„i(to) = —1 and, as it follows [79] 
from (63), P' = cos^ 9. Correspondingly, the average probability takes the form: 

P{ue Ve, ts, to) = 1 - ^ sin^ 261 , (65) 

which is the average two-neutrino vacuum oscillation probability. Thus, if the so- 
lar neutrino transitions are extremely nonadiabatic, the Vg undergo oscillations in 

An expression for the "jump" probability corresponding to the case of density {Ne{t)) 
varying linearly along the neutrino path was derived a long time ago by Landau and 
Zcnor [87]. An analytic description of the solar neutrino transitions based on the 
linear approximation for the change of Ne in the Sun and on the Landau-Zener 
result was proposed in [88]. The drawbacks of this description, which is less accurate 
[83] than the description based on the results obtained in the exponential density 
approximation, were discussed, e.g., in [82], [83], [86]. 
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the Sun as in vacuum. We get the same result, eq. (65), if Ne{to){l — ta,n29)~^ < 
N^^^, i.e., when E/Aw? is sufBciently small so that the resonance density exceeds 
the density in the point of neutrino production. In this case [83] the Ve transi- 
tions are adiabatic (P' = 0) and again the ^ ^fj,(r) oscillations take place in 
the Sun as in vacuum: cos 29m{ta) = cos 29 and P{i'e ^ t^e',ts,to) — 1~ 5 ^in^ 29. 

Let us note that the general aspects of the discussion and the results presented 
above are valid also in the case of solar neutrino transitions into sterile neutrino, 
Ue Vs. In particular, the average probability ^©(fe Ve;ts,to) in this case 
is given effectively by (62) and (63) with [76] Ne{tQ) replaced by {Ne{to) — 
l/2Nn{tf))) in the expression for cos29m{to), Nn{to) being the neutron number 
density of neutrons in the point of neutrino production in the Sun. 

The probability P{i'e Ve;ts,to) is shown as function of E/Am^ for three 
values of sin^ 26 = 0.8; 0.2; 5 x 10"^ in Figs. 7a - 7c. 

Further details concerning the analytic description of the matter -enhanced 
transitions of the solar neutrinos in the Sun can be found in [79], [82]-[86], [88], 
[89]. Exact analytic results for the probability of various possible two- neutrino 
matter-enhanced transitions in a medium {vp — !■ v^i,(t) or the inverse, Vp ^fi,(f) 
or the inverse, v^. ^ij,{t) or the inverse, i/g, etc.), which are based solely 

on the general properties of the system of evolution equations (32) (and do not 
make use of the explicit form of the functions e{t) and e'{t)) are given in [89]. 

Earlier studies (from 1993 - 1994) of the possibility to explain the solar 
neutrino problem in terms of the hypothesis of matter-enhanced Uf. — > ^'^^(t) 
transitions of solar neutrinos have shown [59] that the data admits, in general, 
two types of MSW solutions: a small mixing angle nonadiabatic solution for 
10~^ < sin^ 29 ^ 10~^, and a large mixing angle adiabatic one for approximately 
0.60 ^ sin^ 29 <, 0.95, with the allowed values of Am^ lying in the interval 
10""^ eV^ ^ Am'^ <, lO""* eV^. The terms "nonadiabatic" and "adiabatic" refer 
to the type of transitions the neutrinos undergo in the corresponding cases. 
It was also shown (sec, e.g., [56], [61]) that in the case of Vs transitions 

only a small mixing angle nonadiabatic solution, analogous to the ^ij.(t) 
nonadiabatic solution, is allowed by the data. 

Recently the MSW solutions of the solar neutrino problem have been re- 
examined [72], [68], [69] (exploiting the x^— method) using the data (1), (2), 
and (4), the GALLEX result from 51 runs of measurements, RGALLEx(Ge) = 
(69.7 ± 6.7t 4 5) SNU, and the Super-Kamiokande result from 201.6 days of 
measurements (~ 3000 events), #1^ = (2.65 to^^^ to.io) x 10*^ Ve/cra^/sec . 
The analysis was based on the predictions of the solar model of [29] with heavy 
element diffusion for the electron and neutron number density distributions" and 
for the relevant pp, pep, ''Be, ^B and CNO components of the solar neutrino flux. 
The uncertainties in the predictions for the fluxes estimated in [29] as well as the 
uncertainties of the different solar neutrino detection reaction cross-sections were 



All solar models compatible with the currently existing observational constraints (he- 
lioseismological and other) predict practically the same electron and neutron number 
density distributions in the Sun. 
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Fig. 7. (c)The solax i>e survival probability [90]Po(!Ve —^ i^e;is,to), Eq. (62), averaged 
over the region of production the pp (solid line), pep (long-dash-dotted line), '^^N 
(dashed line), ^Be (dash-dotted line), ^^O (long-dashcd line) and (dotted line) 
neutrinos for sin^ 261 = 0.8 (o); 0.2 (6); 0.005 (c) as a function of E/Am^. Figures a 
and b correspond to Ve — > ^'/^(r) transitions, while figure c corresponds to Ve — ► fa 
transitions. 
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taken into account. The probability pQive — > ^e',ts,to) was calculated following 

the prescriptions given in [83] . The results obtained in the cases of Vg — > ^/^(t) 
and — > i>s transitions are depicted in Figs. 8 and 9, respectively. 

The solid line contours in Fig. 8 denote regions allowed by the data from 
the Homestake, Kamiokande, SAGE and GALLEX experiments, while the dark 
shaded regions have been obtained by including the Super-Kamiokande data in 
the analysis. Thus, the dark shaded areas represent the regions allowed by the 
mean event rate data from all experiments. The solid line contours in Fig. 9 
denote the region allowed by the data (at 95% and at 99% C.L.). 

The current solar neutrino data are best described assuming the solar neu- 
trinos undergo small mixing angle i>e ^fj.{T) matter-enhanced transitions [72], 
[68], [69] (for this nonadiabatic solution one has Xmin = 0-^ d.f.)). The quality 
of the fit of the data is somewhat worse in the case of the large mixing angle or 
adiabatic solution {Xmin somewhat larger: Xmin = 1-5) ■ A similar quality of 
the fit of the data is provided also by the hypothesis of transitions into a sterile 
neutrino, i^e t^s, at small mixing angles (Xmin — ^■^)- contrast, the large 
mixing angle i^^ >^s transition solution is practically excluded as a possible 
explanation of the solar neutrino deficit [72], [68], [69] (it is ruled out at 99.98% 
C.L. (xL„ = 21, 3 d.f.) by the data). 

The values of the parameters Am^ and sin^ 29 for which one obtains (at 95% 
C.L.) the small mixing angle i^e '^/j-ir) transition solution of the solar neutrino 
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Fig. 8. Regions of values of tlie parameters Arri^ and sin'^ 29 {the black areas) for which 
the matter-enhanced — > ^'^{t) transitions of solar allows to describe (at 95% C.L.) 
the solar neutrino data (from [72]). For further details see the text. 



problem lie in the region: 



3.8 X lO^^eV^ <, Am^ <, 10"^eV^ 
3.5 X 10"^ < sin^ 29 < 1.4 x 10"^ 



(66a) 
(66&) 

As Figs. 7 and 8 show, the small mixing angle i^e — > solution region is very 
similar in shape and magnitude to the region of the Ve ^fj-ir) solution, (66a) 
and (66b), but is shifted with respect to the latter by a factor of ~ 1.3 to smaller 
values of Am? . 

We have seen that there can be large uncertainties in the solar model pre- 
dictions for the total flux of neutrinos and that the predictions for the ^Be 
neutrino flux vary by ~25%. The question of how stable are the MSW solutions 
of the solar neutrino problem discussed above with respect to changes in the pre- 
dictions for the two fluxes and ^bc naturally arises. A rather comprehensive 
answer to this question for the ^^^(T) transition solution was given in [60], 

and for the solution with v,, transitions into a sterile neutrino, Ve Vs - in [59]. 
These studies showed, in particular, that the existence of the MSW solutions of 
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Fig. 9. Allowed region of values of the parameters Am^ and sin'^ 29 corresponding to 
the matter-enhanced Ve —> i^s transition solution of the solar neutrino problem (from 
[72]). For further details see the text. 

the solar neutrino problem is remarkably stable with respect to variations in the 
predictions for the *B and ''Be neutrino fluxes. 

6.3 A Detour: MSW Transitions of Solar Neutrinos in the Sun and 
the Hydrogen Atom 

As we have indicated, the two-neutrino matter-enhanced i^^ '^m(t) transitions 
of solar neutrinos at small mixing angles provide the best description of the 
solar neutrino data. In the present subsection we demonstrate [89] that the 
second order differential equation for the probability amplitude Af^{t, tg) of solar 
r^e survival coincides in form in the case of solar electron number density iVg (i) 
changing exponentially along the neutrino path, Eq. (50), with the Schrodinger 
equation for the radial part of the non-relativistic wave function of the hydrogen 
atom, and we comment briefly on this interesting coincidence. 

Using the first equation in (32) to express A^(i, to) in terms oi Ae{t, to) and its 
time derivative, which gives Afj,{t,to) = ^ {e{t) + i-^)Ae{t,to), and substituting 
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■^ij.{t,to) thus found in the second equation in (32), we obtain a second order 
differential equation for Ae{t,tQ): 

I ^ + [e' + e"-ie]|^e(t,to) = 0, (67) 

where e = and e{t) and e' are given by (33). Introducing the dimensionless 
variable 

Z = iro\/2GFNe{to)e-'-^, Zo = Z{t = to), (68) 
and making the substitution 

A,s,to) ^ A{v, I.,) - iz/z„r- ,-(^-^o)+i/;e(t')<it' ^gg^ 

we find that the amplitude A'^{t,to) satisfies [79], [80], [84] the confluent hyper- 
geometric equation [81]: 

where [84] 

Am"^ .2/1 1 . Am'^ . 

a = l + tro -^sm e, c=l + tro-^- (71) 

Equation (70) coincides in form with the Schrodingor (energy eigenvalue) 
equation obeyed by the radial part, tpkiir), of the non-relativistic wave func- 
tion of the hydrogen atom [91], >f'(r) = ^■ijjki{r)Yi„i{9' ,(p'), where r, 6' and (i' 
are the spherical coordinates of the electron in the proton's rest frame, / and 
m are the orbital momentum quantum numbers (m = —I, I), k is the quan- 
tum number labeling (together with 1) the electron energy," Eki {Eki < 0), 
and Yim{0' are the spherical harmonics. To be more precise, the function 
'tpl-liZ) = Z~'^l'^ e^/^ '^klif) satisfies equation (70), where the variable Z and 
the parameters a and c are in this case related to the physical quantities char- 
acterizing the hydrogen atom: 

Z = 2—^/-Eki/Ei, a = aki=l + l- V-Ei/Eki, c=ci = 2{l + l). (72) 
do 

Here oq = %/{mf.e'^) is the Bohr radius and Ej = mee^/{2h'^) = 13.6 eV is the 
ionization energy of the hydrogen atom. It is remarkable that the behaviour of 
such different physical systems as solar neutrinos undergoing matter-enhanced 
transitions in the Sun and the non-relativistic hydrogen atom are governed by 
one and the same differential equation. 

The properties of the linearly independent solutions of equation (70), i.e., of 
the confluent hypergeometric functions, 'P{a,c;Z), as well as their asymptotic 
series expansions, arc well-known [81]. Any solution of (70) can be expressed 
as a linear combination of two linearly independent solutions of (70), ^{a, c; Z) 

° The principal quantum number is equal to {k + I) [91]. 
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and (?(a — c + 1,2 — c; Z), which are distinguished from other sets of hn- 
early independent confluent hypergeometric functions by their behaviour when 
Z ^ 0: ^{a',c';Z = 0) = 1, a',c' ^ 0,-1,-2,..., a' and c' being arbitrary 
parameters. Explicit expressions for the probability amplitudes A{h',. v^) and 
A{ve '^ix{t)) ill terms of the functions ${a,c;Z) and <P{a — c + 1,2 — c; Z) 
were derived in [84], [92]. In the case of MSW transitions of solar neutrinos 
{Ne{ts) = 0) these expressions have an especially simple form: they arc given by 
the corresponding vacuum oscillation probability amplitudes "distorted" by the 
values of the functions ??(a', c'; Z) in the initial point of the neutrino trajectory, 

^K^i^^(r)) = ^ sm29 {^(a-c,2-c;Zo)-e*(*-*°)^ ^(a-l,c;Zo) }, 

(73) 

etc., where Zq, a and c are defined in (68) and (71). In the limit \Zo\ — > 0, which 
corresponds to zero electron number density, expression (73) reduces (up to an 
irrelevant common phase factor) to the one for oscillations in vacuum, Eq. (23). 

It is well-known that the requirement of a correct asymptotic behaviour of 
the wave function V-'fei(^) at large r leads to the quantization condition for the 
energy of the electron, Eki, in the hydrogen atom [91] : Eki = —Ei/{k + l)"^, 
{k+l) = 1, 2, ... (/ = 0, 1, 2, {k+l) — l). Technically, the condition is derived by 
using the asymptotic scries expansion of the confluent hypergeometric functions 
in inverse powers of the argument Z [81] (one has Z ^ oo when r ^ oo, see (72)). 
The same asymptotic series expansion in the case of the solutions describing the 
MSW transitions of solar neutrinos in the Sun (we have \Zq\ ^ 520 in this 
case [84]) permitted to derive i) the simple expression for the relevant "jump" 
probability [79] P', Eq. (63), and ii) explicit expressions for the oscillating terms 
in the solar v^. survival probability [84]. Expression (63) is a basic ingredient of 
the most precise simple analytic description of the two- neutrino matter-enhanced 
transitions of solar neutrinos in the Sun, available at present [83]. 

7 The Solar Neutrino Problem: Outlook 

After being with us for ^25 years the solar neutrino problem still remains un- 
solved. With the accumulation of the quantitatively new data provided by the 
Ga-Ge experiments the problem acquired a novel aspect: the constraints on the 
^Be neutrino flux following from the data imply a significantly smaller value of 
than is predicted by the solar models. The data of both Davis et al. and 
Kamiokande experiments have to be incorrect in order for the indicated conclu- 
sion to be not valid. The vacuum oscillations and MSW transitions of the solar 
neutrinos continue to be viable and very attractive solutions of the problem. 

The start of the Super-Kamiokande experiment on April 1, 1996, and the 
presentation of the first preliminary data from this experiment at the "Neutrino 
'96" International Conference in June of the same year [22], marked the begin- 
ning of a new era in the experimental studies of solar neutrinos. This is the era 
of high statistics experiments with real time event detection and capabilities to 
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perforin high precision spectrum, seasonal variation [2], [56], day-night asym- 
metry (see, e.g., [93], [94] and the articles quoted therein), etc., measurements. 
Such capabilities are of crucial importance, in particular, for understanding the 
true cause of the solar neutrino deficit. 

The preceding period 1967 - 1996 of solar neutrino measurements, which 
began when the epic Homestake (Cl-Ar) experiment started to collect data 
[15], [17], is marked by several remarkable achievements which, given their scale 
and the time and the efforts they took, make this period rather an epoch. For the 
first time neutrinos emitted by the Sun have been observed. The thermo-nuclear 
reaction theory of solar energy generation was confirmed by the detection by 
GALLEX and SAGE experiments of the lower energy solar neutrinos produced 
in the corresponding fusion nuclear reactions. More generally, this result con- 
firms a fundamental aspect of the theory of stellar evolution regarding the role 
played by the nuclear fusion reactions. Finally, the solar neutrino data gathered 
in the indicated period provided, when compared with the predictions of the so- 
lar models, indirect evidences for an "unconventional" behaviour (e.g., vacuum 
oscillations, and/or matter-enhanced transitions, etc.) of the solar neutrinos on 
their way to the Earth. This in turn is the strongest indication we presently have 
for the existence of new physics beyond that predicted by the standard theory 
of electrowcak and strong interactions. 

The Super-Kamiokande is the first operating of a group of new generation 
detectors, SNO [95], BOREXINO [96], ICARUS [97], HELLAZ [98], etc., which 
will allow one to perform more d(rt,ailcd and accurate studies of the solar neutrino 
flux reaching the Earth. As is well known, Super-Kamiokande, SNO and ICARUS 
experiments will study the component of the solar neutrino flux at energies of 
solar neutrinos E ^ (5—6) MeV; the BOREXINO detector is designed to provide 
information about the 0.862 MeV ^Be component of the flux: approximately 90% 
of the signal produced by the solar neutrinos in the BOREXINO detector (~50 
events/day according to the reference model [29]) is predicted to be due to the 
^Bc- neutrinos. The HELLAZ apparatus is envisaged to measure the total flux 
and the spectrum of the pp neutrinos ^ in the energy interval E = (0.22 — 
0.41) MeV. 

The SNO experiment is expected to begin to take data in 1998. The construc- 
tion of the BOREXINO detector is under way and is planned to be completed 
by the end of 1998. A prototype of the ICARUS apparatus has been successfully 
tested and the construction of the first 600 ton module has started. The feasi- 
bility studies for the HELLAZ detector have been intensified with the building 
of a small prototype at College de France [98] . Our aspirations to find the cause 
of the solar neutrino deficit established by the results of the spectacular solar 
neutrino experiments of the first generation [15], [17], [19]-[21], and confirmed 
by the first results from the Super-Kamiokande detector, and to get additional 
independent information about the physical conditions in the central part of the 
Sun, are presently associated with the more precise and diverse data the second 



The HELLAZ detector can be utilized for studies of the ^Be neutrino flux as well. 
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generation detectors are expected to provide. All these are planned to be high 
statistics (typically ^^3000 solar neutrino events/year, Super-Kamiokande is ex- 
pected to collect ^10000 events/year), i.e., high precision, experiments with real 
time event detection. 

In SNO experiment the *B neutrinos will be detected via the charged current 
and the neutral current reactions on deuterium: fe + D — * e~ + p + p, and v + 
D — > + p + n; the measurement of the kinetic energy of the electron in the first 
reaction will permit to search for possible deformations of the spectrum of *B 
neutrinos at S > 6.44 MeV, predicted to exist (sec, e.g., the first article quoted 
in [59] as well as [56], [61]) if solar neutrinos take part in oscillations in vacuum 
on the way to the Earth and/or undergo matter-enhanced transitions in the Sun. 
High precision searches for spectrum deformations will be performed also in the 
Super- Kamiokande experiment in which the energy of the recoil electron from 
the V — Q~ elastic scattering reaction will be measured with a high accuracy. 

The high statistics these experiments will accumulate, the measurement of 
the spectra of final state electrons with the SNO and Super Kamiokande detec- 
tors, and of the ratio of the charged current and the neutral current reaction 
rates with the SNO detector, will make it possible to perform various critical 
tests (see, e.g., [56], [61], [93], [94]) of the vacuum oscillation and the MSW, 
as well as of the other possible neutrino physics solutions [62]-[64], [66] of the 
solar neutrino problem. We may be at the dawn of a major breakthrough in the 
studies of solar neutrinos. It is not excluded, however, that the data from the 
BOREXINO and HELLAZ detectors may be required to get an unambiguous 
answer concerning the cause of the solar neutrino problem [60], [57], [61]. 
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